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1.  INTRODUCTION  AND  SUMMARY 

An  enduring  problem  in  electromagnetic  theory  is  that  of  calculating 
the  response  of  a  horizontal  wire  over  earth  to  incident  electromagnetic 
fields.  An  approximate  solution  of  the  problem  that  is  widely  used  is 
obtained  by  regarding  the  wire  over  earth  as  a  transmission  line. 

The  basic  unknown  of  the  transmission-line  approach  is  the  approxi¬ 
mate  relationship  (from  field  theory)  between  the  line  waves  and  the 
electromagnetic  fields  in  the  earth;  this  relationship  enters  the 
transmission-line  formulation  as  a  surface  impedance  for  the  ground 
plane. 

The  purpose  of  this  paper  is  to  present  a  heretofore  overlooked 
model  formulation  and  exact  solution  to  the  wire-over-ground  surface- 
impedance  problem.  The  surface-impedance  model  of  section  2  results 
from  the  joining  or  transmission-line  waves  to  purely  TM  (transverse 
magnetic)  electromagnetic  soil  fields.  More  precisely,  we  also  show 
(sect.  4)  that  the  TM  surface  impedance  arises  from  a  purely  vector- 
potential  formulation  of  the  wire-over-ground  problem,  if  the  model 
assumed  in  section  2  is  granted.  It  is  then  possible  to  derive  (1)  the 
line-inductance  parameter,  (2)  the  "charge"  distribution,  (3)  the  finite 
width  of  dominant  wave  localization  over  the  ground  plane  (derived 
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consistently  with  the  line-capacitarice  parameter),  and  (4)  the  surface- 
impedance  boundary  relations,  all  from  the  vector-potential  solution. 
The  zero  possibility  of  certain  TE  (transverse  electric)  field  compo¬ 
nents  (i.e.,  all  omitted  field  components)  follows  automatically  from 
the  vector-potential  solution.  The  vector-potential  model  is  indistin¬ 
guishable  from  the  commonly  accepted  definition  of  a  surface-impedance 
loaded  transmission  line.  This  paper  does  not  consider  the  extent 
to  which  higher  order  waves ,  excluded  from  the  TM  model ,  compromise 
the  TM  line  response  (predominantly  at  high  frequencies,  where  the 
transmission-line  "radiation  resistance"  is  large). 

Our  surface  impedance  does  not  agree  with  either  the  low-frequency 
surface  impedance  in  the  work  of  Carson1  or  the  essentially  identical 
low-frequency  surface  impedance  that  is  a  limiting  case  of  the  Wait2 
formulation  of  the  wire-over-ground  problem.  At  higher  frequencies  we 
find  a  guided  surface  wave  ar.d  surface  impedance,  unlike  the  results  in 
the  paper  of  Kikuchi.^ 

Section  3  is  devoted  to  the  two-dimensional  inhomogeneous  plane 
waves  resulting  in  the  surface-impedance  model.  The  connection  between 
the  waves  and  the  complex  refraction  angle  formalism  is  discussed. 


1 J.  R.  Carson,  Bell  Sys.  Tech.  J.,  5_  (October  1926),  539-554. 

2 James  R.  Wait,  Radio  Science,  7_,  No.  6  (June  1972),  675-679. 
.  Kikuchi,  Electrotech.  J.  Japan,  2_,  3/4  (1956),  7  3-78. 
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2.  CALCULATION  OF  Zg,  PER-UNIT-LENGTH  SURFACE  IMPEDANCE 

I  first  describe  the  situation  at  hand  and  give  the  central  re¬ 
sults.  The  derivation  and  discussion  of  the  surface-impedance  model  is 
then  taken  up. 

A  horizontal  wire  of  radius  a,  lying  at  height  h  over  planar  earth, 
is  regarded  as  a  transmission  line  in  essentially  the  same  sense  as  the 
Beverage  or  wave  antenna.1*  The  wire  is  highly  conducting  and  character¬ 
ized  by  a  prescribed  ( ps.r-un it- length )  internal  impedance,  Z^,  where 

=  R^  +  iwJ the  resistance,  R^ ,  and  inductance,  are  functions  of 
radian  frequency  (w  =  2irf),  and  f  is  the  usual  cw  frequency  parameter  in 
hertz . 

The  earth  is  characterized  by  a  scalar,  frequency-dependent  complex 
dielectric-response  function,"’  e(w)  where  e(ui)  =  e^io)  -  ie2(w).  The 
earth  conductivity  is  a(u)  =  ooe2(u>)  =  °g  +  we2r^a)^'  w^*ere  °q  is  3  con_ 
stant  dc  conductivity  and  e  is  the  remaining  loss  part  of  z^. 

iu)t 

I  take  here  an  e  Fourier  convention;  therefore ,  the  line  propaga- 
tion  constant®  is  y  =  a  +  i3,  a(u>)  is  the  attenuation  part,  and  v(u>)  = 
w/BCw)  relates  6  to  the  phase  velocity,  v.  The  square  of  y  is 

^See,  for  example/  J.  A.  Schelkunoff  and  H.  T.  Friis,  Antennas 
(Theory  and  Practice) ,  John  Wiley  and  Sonsr  Inc.  (1952). 

5J .  D.  Jackson,  Classical  Electrodynamics ,  second  edition,  John  Wiley 
and  Sons,  Inc.  (1975),  p  309. 

®R.  W.  P.  King,  Transmission  Line  Theory,  Dover  Publications,  Inc . 
(1965),  p  91 . 
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(1) 


=  iwcfz  +  Z.  +  ibil  )  , 

'  g  i  e ' 

where  l is  per-unit-length  external  inductance  and  C  is  per-unit-length 
capacitance.  The  values  of  C  and  are  exactly  the  values  calculated 
for  a  wire  over  a  perfectly  conducting  ground  plane  (sect.  4).  I  assume 
no  shunt  conductance  (G  =  0)  for  the  wire. 

The  wire  as  a  transmission  line  is  complete  once  we  know  the  per- 
unit-length  surface  impedance  for  the  earth  presence,  namely 
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$2  =  tan'1  (F./Fj  , 


(7) 


w  =  2nh  .  (10) 

Thus,  if  we  substitute  2  of  equation  (3)  into  (1)  and  extract  the 
positive  root6  y  =  a  +  ifi,  we  completely  define  the  uniform 
transmission-line  parameters  with  respect  to  the  basis  waves  from  which 
the  well-known  response  solutions^'®  for  the  terminated,  finite-length 
line  with  arbitrary  excitation  are  constructed. 

The  transmission-line  voltage  is  V  with  a  current  I  on  the  upper 
wire.  It  is  helpful  to  regard  the  line  response  as  a  guided  surface 

6R.  W.  P.  King,  Transmission  Line  Theory,  Dover  Publications,  Inc. 
(1965),  p  91 . 

7P.  F.  Gray,  Nuclear  Electromagnetic  Pulse  Simulation  by  Point  Source 
Injection  Techniques  for  Shielded  and  Unshielded  Penetrations ,  Harry 
Diamond  Laboratories ,  HDL-T R-1737  (December  1 975), 

8See,  for  example,  A.  A.  Smith,  Jr.,  Coupling  of  External  Electro¬ 
magnetic  Fields  to  Transmission  Lines,  Wiley  Interscience  (1977). 
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wave  and  to  first  focus  attention  on  the  earth-surface  magnetic  field. 
For  brevity  and  compactness  of  expression  I  have  decided,  in  most  of 
what  follows,  to  write  down  only  the  frequency-domain  representatives  of 
responses  (denoted  by  underlined  symbols;  all  responses  not  underlined 
denote  real  responses)  and  to  suppress  the  explicit  conversion  of  the 
underlined  quantities  to  real  responses  (unless  confusion  would 
arise).  In  our  right-hand  coordinate  system,  x  is  upward  from  the 
ground  and  z  is  along  the  transmission  line.  For  calculating  per-unit- 
length  Z in  the  frequency  domain,  it  is  convenient  to  suppose  an 
unbounded  line  with  a  Dirac  impulse  point-source  line-voltage  generator 
at  z  =  0  and  to  single  out  the  +z-directed  Fourier  transmission-line 
traveling  waves. 

The  TM  2^  follows  from  a  single  model  assumption  of  simultaneous 
contributions  from  all  y  locations  of  surface  magnetic  field  to  the 
transmission-line  return  current  -I,  namely  at  x  =  0,  and  some  fixed  z, 

-I(t)  =  f  dy  H  ( 0 ,y ,t )  . 

*  «-co  y 

Hence,  in  general,  the  y  integral  of  must  exhibit  common  time  depend¬ 
ence  or 

Hy(0,y,t)  =  h1(t)h2(y)  , 
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r°° 

where  h,  and  are  two  real  functions.  It  follows  from  I  dy  h_ 

12  J  -a>  J  2 

immediately  that  space  frequency  k^  is  zero  in  the  line  waves.  We  may 
then  restrict  ourselves  again  to  h^(u,z)  in 


-I  =  fL  dy  V0,Y) 


wH  exp(-Yz) 
-yg 


(ID 


where  H  is  a  frequency-dependent  constant  and  width  parameter 
,00 

w  =  J  ^  h  (y)  dy. 


The  soil  fields  in  the  limit  x  =  0-  must  also  have  k^  =  0  and  share 
common  h2(y)  and  h^  since  Hy(0,y)  is  the  prescribed  soil  field  boundary 
Hy  distributed  at  x  =  0  over  all  y  instantaneously,  and  all  soil  fields 
are  derivable  from  Hy.  In  (11)  we  could  prescribe  heuristically  from 
magnetostatics  the  results 


where 


and 


Hi  = 


H 

-yg 


H 

-yg 


H  (0,0) 

Zl _ 


2 


h2<y> 


2hw 

2n(h2  +  y2) 


,O0 

hence,  J  h  (y)  dy  =  w  and  h  (0)  =  2w/27ih  =  2  from  equation  (10).  In 
section  4  we  show  that  the  prescription  is  not  necessary  since,  if  the 
prescription  is  a  bit  further  elaborated,  it  has  the  same  results  as  the 
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vector-potential  calculations*  In  the  earth  only  the  TM  electromagnetic 


fields  E  ,  E  ,  and  H  are  kept,  whence 
x  z  y 


E 

-x 


i we  3z  ' 


(  12) 


E 

-z 


iu)£  3x  ' 


H  = 

-y 


ia)yQ  \3x 


(13) 


(  14) 


-►  -► 

•ik*x 

On  taking  the  exponential  dependence  e  ,  if  we  substitute  equations 

(12)  and  (13)  into  (14),  we  obtain 

2  2  2 

u>  pneH  =  fk  +  k  Jh 
0  -y  v  x  z'-y 


or 


where  ky  =  0 . 


+2  2 
k  =  p  Qeoj 


Since  e  is  a  scalar,  k»E  =  0  is  the  statement  of  div  D  =  0  in  the 
soil.  Normal  D  must  be  continuous  at  x  =  0;  thus,  the  connection  at  x  = 
0  of  line  wave  to  soil  waves  involves  only  the  integrated  (defined) 

r  OO  -f  A 

charge  per  unit  length  g  =  dy  D*x  taken  with  respect  to  upward 

normal  x.  Or,  if  we  use  equation  (12)  to  evaluate  g,  we  obtain 


g(x  =  0 ) 


00 

—00 


1 


iu> 


(15) 
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which  is  just  the  continuity  equation 


0 

iwq(  x  =  0)  =  -  —  (-1)  =  -itoCV  ,  (15a) 

3  3z  -  - 

where  equation  (11)  has  been  used  in  (15)  (the  right  side  is  also  the 
time  derivative  of  eq  (72)).  (Remember  that  (j^  g  =  -g(x  =  0).) 

We  need  now  only  superpose  Fourier  waves  in 
+  h2(y)  °° 

-y(x)  =  lir —  dkz  *0^kz^exp^_ikxx^  exP("ikzz)  '  M6) 

where 

,  /  2  2-|1/2  .  , _ , 

kx  =  "(V®  "  V  '  (17) 

+2  2 

constrained  by  k  =  UgGW  ;  the  square  root  in  equation  (17)  taken  in 
quadrant  4  yields  downward-traveling  damped  waves  in  the  earth. 

We  first  place  x  =  0  in  equation  (16),  integrate  (16)  over  dy  from 

ik'z 

-»  to  °°,  and  then  equate  it  to  (11);  we  then  multiply  by  e  and 

integrate  over  all  z,  using  the  following  representation  of  the  Dirac 
delta  function  <5, 

-ifk  -k'lz 

1  /2tt  f*  dz  e  2  z  =  6(k  -  k') 

J  — OO  V  Z  Z ' 


All  this  produces 
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Therefore,  also,  wZ  =  E  (0,0)/H  (0,0)  =  -k  /ioe.  The  solution 

cf  -2  -y  x  9 

follows  from  inserting  equation  (1)  into  (21)  squared.  The  details  are 
given  in  appendix  A. 

We  must  clarify  an  important  point  concerning  the  total  (reaction) 
electric  field  at  x  =  0.  A  derivation  of  the  transmission-line  voltage 
equation  for  the  case  at  hand  gives 

3V  . 

=  -ifz  +  Z.  +  imi  1  +  ElnC(h,z)  -  EinC(0,z)  ,  (22) 

dz  —  g  l  e  —  z  —  z 

which  includes,  for  generality,  a  distributed  electric-field  excitation 
(the  wire  is  assumed  perfectly  thin).  Equation  (22)  follows  from  the 
well-known 

§(E  +  ElnC)*d t  =  0  (23) 

taken  about  a  rectangular  path,  of  thickness  Az,  over  x  =  0  to  h.  We 
first  ignore  and  proceed  with  the  usual  Taylor  expansion  by  taking 
from  (23) 

§E»  d£  =  -§E'i'nC,dt  =  -  "applied  voltage"  .  (24) 

For  each  path  segment,  we  also  define  a  reaction  voltage  for  path  a  -*•  b 

as  V  =  -fk  E*d£.  Thus,  by  invoking  only  the  sign  of  d$,  we  correctly 
ab  'a 

add  to  equation  (24)  as  yet  unspecified  electric  fields  along  the  paths, 
inc 

We  ignore  E^  and  invoke  Z^I  =  Ez(h,z)  and  Z  (-1)  =  Ez(0,z),  according 
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to  my  current  sense  convention  on  choosing  a  clockwise  traverse  of  the 
paths.  Finally,  we  add  to  the  right-hand  side  of  (22)  (shown  completed) 
the  contribution 


3V' 

-iw*  I(Az)  =  -§E'*d£  =  -f-  (Az)  ,  (25) 

e-  —  a  z 

where  V  is  the  voltage  and  E  is  the  electric  field  contributions  in 

linear  superposition  due  solely  to  the  external  inductance  The 

important  point  is  that  the  expression 

E'  =  0  (26) 

z 

in  the  horizontal  paths  of  (25)  is  the  correct  approximation  in  general, 
regardless  of  and  Zg  being  nonzero.  For  this  reason,  E^  in  equation 
(21)  is  the  correct  superposition  component  to  total  reaction  E^ 
(x  =  0). 


To  show  explicitly  that  equation  (26)  is  correct,  we  need  only  exam¬ 
ine  the  signless  quantity 


I N'  |  =  Az  /*;  dt  S'*  — 
1  0  ,  ->■ 


S'l 


9e' 


/rAz  dz  \\  dt  Jl  dy  /;  dx  3/  Hy  ( 0  ,y ) 


(27) 


for  an  infinitesimally  thick  transverse-line  slice,  where  S'  is  a 
Poynting  flux  contribution  due  solely  to  the  E'  of  (25).  Using  the 
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earlier  voltage  definition,  substituting  equations  (11)  and  (25)  into 


(27),  and  integrating  by  parts  with  I(O-)  =  0,  we  find 


l  I  ( t ) Az 


IN'  |  = 


this  shows  that  |N"VAz|  is  indeed  the  correct  magnetic  energy  per  unit 
length  for  the  line  in  the  presence  of  full  conductor  impedance.  Of 
course,  we  redefine  t  =  0-  for  any  plane  selected,  but  the  result  is 
true  in  complete  generality  for  the  arbitrarily  excited  and  terminated 
line,  since  total  S'  is  linear  in  the  superposition  voltage  component 
V'. 


Because  of  equation  (25)  it  is  clear  that  in  a  consistent  approxima¬ 
tion,  the  product  J^w  must  act  instantaneously;  hence,  both  and  w 
must  be  real  and  constant  (independent  of  frequency).  Parenthetically, 
C  (through  the  upper  medium  dielectric  constant  e  )  can  depend  on 

frequency  {i  C  =  y  e  (uj),  where  y  is  a  real  constant),  but,  for  an 
e  up  up  up 

air  medium,  the  effect  of  e^fto)  is  insignificant. 

Qualitatively,  the  earth  in  this  model  is  infinite  and  spatially 
homogeneous,  so  that  no  upward-going  waves  originate  in  the  soil  vol¬ 
ume.  In  the  air  the  dominant  waves  are  the  one-dimensional  line  waves 
with  k  =  k^z  .  Higher  order  waves  bring  in  k^,  wave-vector  components 
perpendicular  to  k^z. 
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We  can  roughly  estimate  when  these  higher  order  waves  containing  k^ 
in  the  air  will  be  important.  To  do  this  we  can  borrow  an  argument  from 
Kompaneyets9  which  originally  applied  to  limits  of  accuracy  in  radio¬ 
location.  The  argument  is  that  one  attempts  to  confine  these  higher 

order  waves  roughly  in  a  half  cone  flared  about  k  z,  so  that  the  waves 

z 

contribute  to  a  forward  directivity  of  the  "antenna."  Similar  isosceles 
triangles  are  applicable  with  height  or  B  and  half  base  length 

Wg/2  or  k^;  hence,  the  wave  confinement  is  estimated  at 
Wg/2£g  =  k^/3,  and  the  space-frequency  wave-uncertainty  relation 
is  4(wg/2)k^  >  2ir.  Eliminating  k^  we  find 

(6wq)2  >  4n2(ZQ/X)  , 

where  3  =  2tt/X,  or 

(3h)2  >  rr2(ZQ/x)  ,  (28) 

if  Wq/2  =  h  is  taken  as  the  appropriate  estimate.  For  fixed  l Q  versus 
Iq/X,  Iq/X  =  0.25  is  roughly  the  onset  of  the  higher  order  waves,  and 
Zq/X  =  1  roughly  the  fully  developed  higher  order  waves.  Thus  equation 
(28)  gives  about 

{ 3h) 2  >2+10  , 

which  indicates  roughly  the  buildup  region  of  the  higher  order  waves. 

9A.  S.  Kompaneyets ,  Theoretical  Physics,  Dover  Publications,  Inc. 
(1962),  p  179. 
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2 

Usually,  in  transmission-line  applications  one  attempts  to  have  (Bh)  << 
1.  For  the  infinite  line,  of  course,  is  any  line  section  of 

length  =  X,  which  is  instantaneously  full  of  the  waves  of  interest 
(steady  waves,  now). 

Any  e(to)  that  exhibits  e2^ei  <  1  at  high  frequency  will  yield  a 
reactance  crossover  (X  =  0),  followed  by  increasingly  negative  X  ,  at 
some  high  frequency.  It  is  easy  to  see  directly  from  the  Maxwell  equa¬ 
tions  that  the  impedance  must  be  a  pure  negative  reactance  if  a  =  0  re¬ 
sults  from  soil  with  =  0  (sect.  3).  Hence,  the  reactance  X^  is  also 
negative  when  is  small.  The  Xg  crossover  can  be  calculated 

from  -j  “  =  0,  using  equations  (5)  and  (7).  This  equation  results 

in 

n6  +  n4(2r  -  1)  +  n2(4r  -  5)  +  2r  -  3  =  0  , 

r  =  (^e  +  *i)c/Vi  ' 

where  =  0  and  n  =  e^/e^.  Usually,  r  is  small,  so  we  can  solve 

6  4  2 

nQ  -  nQ  -  hQ5  -3=0  resulting  in  -  1.73  or  e2/e1  ~  1*7.  Solving 
the  r  *  0  equation  perturbatively  (h  =  +  X),  we  find  roughly  X  = 
-r/3hQ  to  order  n  \  The  zeroth  solution  is  accurate  enough  for  most 
purposes,  though  easily  made  exact  numerically. 
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3.  TWO-DIMENSIONAL  INHOMOGENEOUS  PLANE  WAVES 

From  equations  (17)  and  (19),  the  wave  vectors  k  and  kz  are 

kx  =  -(A  -  iB)  ,  (29) 

k  =  (5  -  ia  . 
z 

Thus,  A  and  B  are  positive  definite  for  the  quadrant  4  square  root  in 

equation  (21).  We  can  solve  equations  (17),  (21),  and  (29)  for  A  and  B 

(the  only  practical  way  to  obtain  these  quantities): 

A  =  usw^Rg  +  e2Xg)  ,  (30) 

B  =  u>w(e  R  -  e_x  )  .  (31) 

^  2  q  i  g 

On  the  other  hand,  we  can  take  the  square  root  in  equation  (21) 
directly,  producing 

A  -  iB  =  (p^  +  P^)1^4  [cos  (r/2  +  't'c/'O  "  *  sin  l1*/2  +  <t>c/2)]  (32) 

where 

♦c  = tan_1  [p/(-ip2n]  * 
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with 


P,  =  VJ„U)^e„  -  u)C(r.  +  R  )  =  2AB  , 

1  2  v  1  g; 

P„  =  m<soj2£.  -  wcTx  +  wf  l  +  i.)l  =  A2  -  B2 
201  L  g  ^  e 

The  root  A  -  iB  of  equation  (32)  lies  always  in  octant  7  of  quadrant  4, 
and  A  and  B  satisfy  the  phase  identity 

-it/2  +  j  tan-1  [2AB/(B2  -  A2 )  ]  =  tan_1(-B/A)  , 

which  holds  also  for  the  other  root  (octant  3)  -A  +  iB.  This  root 

follows  from  the  previous  identity  under  A  -A,  B  +  -B;  the  angles  are 
then  reckoned  in  the  opposite  sense  (positive  clockwise  from  the  -x 

axis) . 

The  unique  situation  here  is  that  this  model  excludes  octant  8  of 
quadrant  4  since  the  e2/e1  >>  1  corresponds  to  the  A  =  B  limit 

between  octants  7  and  8,  approaching  from  octant  7  or  A  <  B.  As  0 

the  other  boundary  of  octant  7  is  reached. 

In  the  lossless  limit  (e„  =0,  R.  =0),  one  finds  R  =  0,  a  =  0, 

2  l  g 

A  =  0,  <f>1  =  -it/2,  <j>  =  0,  and  lossless  B  =  Bp,  or 

B0  =  jHXg  +  “(*e  +  Ai)l  “  V/)T/4  '  (33) 
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This  limit  is  correct  because  the  divergence  of  the  real  part  of  the 
complex  Poynting  vector  must  vanish  in  the  absence  of  o  losses;  hence, 
Rg  must  vanish  so  that  a  vanishes  and  a3  +  AB  =  0. 

Clearly,  the  branch  is  unique  since  Zg  is  correct*  for  e  =  e  ^ ,  and 
the  analytic  continuation  of  Zg  to  e  =  -  ie^  is  the  desired  answer. 
In  order  for  AB  to  be  always  positive,  there  is  a  formal  restriction 
that  R^  <  oI^/Eq.  Normally,  the  inequality  is  never  violated. 

The  real  soil  at  high  frequency  is  in  the  region  of  <  and 
the  lossless  solution  just  given  is  merely  an  idealization  of  the  real 
case  with  A  small  but  nonzero.  The  qualitative  point  here  is  that  the 
damping  in  the  soil  tends  to  be  purely  reactive  (pure  being  e2  -  0)  when 
e2/e1  is  small,  yet  the  high-frequency  solution  still  exhibits  a  pro¬ 
nounced  skin  effect  with  B- ^  small.  Thus,  the  real  soil  naturally 
favors  this  octant  7  solution  at  both  high  and  low  frequency. 


It  is  evident  that  the  derivatives  3H  /3x  and  3E  /3x  do  not  enter 

y  z 

explicitly  the  joining  of  the  transmission-line  solution  to  soil 

fields.  Tangential  E  =  -Z  I  is  known  and  is  sufficient  to  predict  the 

-z  g- 

purely  transmitted  fields  into  the  soil.  It  is  of  interest  now  to  see 


*It  is  readily  shown  that  our 
two  possibilities  that  exhibits  X g 


Zg  solution  for  =  0  is  the  one  of 
real  for  unrestricted  . 
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how  the  wave-vector  formalism  relates  to  complex  ref raction-angle 


formalism  and  also  to  surface  waves  of  the  more  familiar  kind, 
need  to  collect  a  few  results. 


First  of  all,  we  can  write  the  complex  wave  vector  as 


k  =  e0m  -  iV 


where  the  normal  m  to  the  constant-phase  planes  and  normal 
constant-amplitude  planes  are 


m  =  -cos  lj^x  +  sin  z  , 

•»  ' 

n  =  -cos  i^x  +  sin  , 

A  A 

where  x  and  z  are  unit  vectors,  and  angles  and  are 
counterclockwise  from  the  -x  axis.  The  real  wave  number  is 

=  Re(k»m) 

We  must  demand  that  -A  +  iB  =  B„m  -  ia„n  and  -ia  +  8  =  B„m 

Ox  Ox  0  z 


We  first 


n  to  the 


positive 


-  iann  ; 
0  z 


hence,  the  following  four  relations  must  all  hold  simultaneously: 


The  true  phase  velocity  is  therefore 


Consider  now  a  complex  refraction  angle 


The  vector  analogous  to 


mis 


-cos  0^x  +  sin  0^z 


The  phase  of  the  refracted  wave  is 


where 


k*x 


Vi 


->• 

•x 


1/2 


1/2  , 

( n 


i  *c ) 


Here,  n  and  <  are  given  by  the  octant  8  square  root 


n  =  [2^ +  Rei^ 


1/2 


[|(|e|Q  -  nEl)]1/2  , 


,  .  .  (2  2 1 1 /2 

where  |e|  =  [  e  +  z^)  .  Hence, 


2 

n 


2 

K  = 


V 


n<  =  z^l/2  , 


where  SJ  is  a  real  positive  definite  quantity,  to  be  determined  momentar¬ 
ily. 
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The  k  =  (-A  +  iB)x 
form  if 

cos  8 

sin  0 

(where  £*  =  e1  +  i£2)» 

The  formal  Snell  law  is 

or 

One  evaluates  the  0 

cos  0^  =  [ ( 
sin  =  [( 


1 


+  (0  -  ia)z  form  of  k  is  identical  to  the 


Vi 


1/2  1/2 

=  (£*fl)  (A  -  iB)/u0  lelojfl  ,  (36 

1/2  1/2 

=  (e*fl)  (0  -  ia)/uQ  I e |(ofi  (37 

1  /2 

since  (£*ft)  =  n  +  ic  and 

n^  +  =  ft ) £ | 


k„n*z  =  k 
1  z 


sin  0,  =  k  /k„  . 

1  z  1 


1  expressions,  equations  (36)  and  (37),  to  find 


1  /2 

nA  +  <B)  -  i(nB  -  <A)]/u0  IeIwSI  ,  (38 

n0  +  tea)  +  i(<0  -  na)]/\i^2  |£|wfi  •  (39 
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The  factors  in  parentheses  are  positive  definite  since  8/a  >  1  and  n/< 
>  1  ( for  e  ^  >  0 ) • 

Note  now  that  the  dispersion  relation  is 

+2  2+2  2 

k  =  k^n^  =  ft  (40 ) 

+2  +2 
since  n^  =  1.  Using  (38)  and  (39)  above,  the  condition  n^  =  1  in  real 

and  imaginary  equality  gives  the  equation  pair 

O  O  0  O  0  0  0 

e  ^ft  ( A  +  8  -  a  -  B  )  +  2e2ft(a6  +  AB)  =  | e  |  u  ft  ,  (41) 

0  0  0  0 

e  ft(a6  +  AB)  =  (e2/2)ft(A  +6  -  a  -  B  )  ,  (42) 

2  2 

if  we  use  n  -  k  and  n<  above  to  reduce.  In  unknowns  2(AB  +  a8)  and 
2  2  2  2 

A  +6~a~B/  the  solution  to  the  equation  pair  is 

2  2  2  2  2 

A  +6  -a  -  B  =  UqEjW  ft  t  (43) 

2 

2(a6  +  AB)  =  U0e2w  ft  •  (44) 

We  see  immediately  that  if  a8  +  AB  (and  e^)  is  positive  definite,  ft  must 
be  the  positive  definite 
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(45) 


n  =  I  A2  +  e2  -  a2  -  b2\/Vq€^2  . 


Hence,  the  real  equality  equation  (43)  of  the  dispersion  relation  solu¬ 


tion  simply  defines 


A2  +  g2  -  a2  -  B2  =  |  A2  +  B2  -  a2  -  B2| 


or  elimination  to  the  identity  (i.e.,  x  =  |x|  is  true  on  eliminating  the 
left-hand  side,  resulting  in  |x|  =  |x|).  A  proof  of  this  absolute 
condition  is  given  in  equation  (60)  below.  Note  that  one  can  write 


equivalently 


+  ->•  * 
n.  =  n  -  in, 
1  r  I 


with  pure  imaginary  inj  and  pure  real  n^ : 


^  A  A  « 

nr  =  [~(nA  +  <B)x  +  (n0  +  k(x)z]/\1q  |e|tuft  , 


rij.  =  [-(<6  -  na)z  - 


(nB  -  <A)x]/yJ  I e | wfi 


*2 

Then  n1 


=  1  becomes 


*2  +2 

nr  ”  ni  =  1  ' 


+  -*■ 

n  «n,  =0  . 

r  I 
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These  equations  are  identical  to  equations  (41)  and  (42)  when  reduced 


(50) 


|§|/|H!  =  ( cip  +  eJ)1/2/u.|e|  , 


and  second,  also  using  a0  +  AB  =  ao®om*n' 

m*n  =  cos  -  i|»2)  =  a[ 6^  -  .  (51) 

2  2  2  2 

On  imposing  0^  -  ■+■  1 0^  -  ot  |  ,  equation  (51)  is  the  same  as  equation 

(44)  of  the  dispersion  relation.  That  cos  ( -  \p2 )  is  even  means  that 
the  power  spectrum  of  the  ohmic  heating  density  is  the  same  for  left¬ 
going  line  waves,  just  as  for  right-going  waves.  The  left-going  waves 
are  obtainable  by  y  -y  (z  <  0);  hence,  m  and  n  are  mirror  symmetric 
across  the  -x  axis  (i^  and  \p^  are  then  negative  and  clockwise  from  the 
-x  axis).  It  is  evident  that  the  wave-creation  origin  can  be  placed  at 
any  z.  The  validity  of  the  result  Zg  with  respect  to  the  terminated 
line  is  due  to  phenomenological  terminations  being  local  (extensionless) 
and  hence  not  changing  Y(w)  on  reflection  of  the  elementary  basis 
waves.  The  terminations  introduce  purely  temporal  delays  (as  well  as 
amplitude  changes)  between  incoming  and  outgoing  basis  waves  in  the 
finite  line  solution. 

It  is  appropriate  now  to  solve  for  0  and  explain  the  apparent  medium 

+2 

til  aspect  of  the  previous  work.  No  matter  what  k  is  (see  app  B)  ,  one 
must  have  the  identity  (using  eq  (50)) 
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(52) 


|e|2/|h|2  =  (um^/lk^l  =  (a^  +  )/to^  I  e 


2„+2,  /  2  .  „2  -i  ,  2  ,  ,  2 


2  -y  -+  ^2  2 

where  |k  |  r  |k*k|.  On  substituting  |k  |  =  y Q | e |  tu  ft  into  equation 
(52),  and  using  the  definition  of  ft  in  (45)  and  also  in  (35),  we  find 
the  ratio  R  on  squaring  to  be 


R  =  (82  -  Uq)2/(Bq  +  a2Q)2  =  uye2te|2/(«2  +  Bq)4  •  (53) 

We  can  also  complete  the  square  using  cos  from  equation  (51), 

resulting  in 


=  (a?  +  -  («We2)(B2  -  a2)2  sec2  ^  .  (54) 


Forming  R  from  equation  (54)  and  equating  it  to  (53)  gives 


(«2  +  32)4  =  y4o)8|e|4[l  +  (e^/lel^)  tan'1  -  *2)]  ;  (55) 


2  . ,  ,  2 


hence,  from  (52)  we  hav"1 


1 1 1 / I H |  =  jf  1/2(u0/|e|] 


1/2 


(56) 


Using  (55)  we  find  explicitly  that 


2_1/“  =  [l  +  (£2/,e|2)  ta°2  (♦i  _ 


(57) 
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The  vector  norms  are  now 


(  2  ^  2 i 1/2  ,,+.2)1/2  n-1/2  f  ,  ,) 

l«0  V  =  l |k|  J  =  n  “l^olel J 


1/2 


(58) 


-►2  ->-2  2 

Ik^l/jkl2  =  o 


-► 

the  ratio  is  1  for  a  homogenous  plane  wave  having  Re(k)  parallel  to 
Im(k)  (for  example,  SI  approaches  1  arbitrarily  closely  as  e  /e.|  tends 
to  infinity). 


From  the  standpoint  of  equation  (56),  the  transmission-line  surface 
wave  transmits  into  an  apparent  medium  with  spatial  anisotropy  due  to 
in  equation  (57),  despite  the  fact  that  the  bulk  medium  is  iso¬ 
tropic  with  respect  to  e.  The  deviation  of  SI  from  unity  is  typically 
such  that  the  right-hand  side  of  (57)  exceeds  unity  by  about  5  percent 
or  so  at  the  higher  frequencies. 


Clearly,  the  bulk  medium  e  is  still  e  and  not  fie.  To  see  this  we 
can  use  o  as  an  example  and  examine  heating  per  unit  volume  or 


f0  J'B  at  =  V2  /0 


®  dll)  .  +  .  2 

27  0,5| 


Using  the  k  norms  above, 


however,  along  with  (48)  we  find 
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|E|2  =  |k  |  2  |h(  2/o>2|  e  |  2  =  (u0/fl|e|]|H 


and  also 


1 12  =  |k2||H|2/|c|Vfi2  =  ( I  k2  |  /uQ  t  e  |w2ft)  ( I  e  I )  |h 


t,2 


with  the  underlined  factor  equal  to  1  in  the  last  equation.  Thus, 


>  2 

substituting  |E|  into  the  integral  exhibits  only  the  partial  apparent 


medium  aspect  of  equation  (56). 


The  partial  quadrature  solution,  entirely  in  terms  of  aQ ,  BQ  ,  ipf , 
and  proceeds  from  the  norm  relation 


a2  +  e2  =  v  |  e  I  u)2Sl~ 1 
0  0  0 


and  the  absolute  value  of  equation  (51),  which  becomes 


O060  =  V“°  S6C  (*1  "  ^2^/2 


solved  by  square  completion.  We  find 


1  2 

a0  +  80  =  20  "  1“  Uo,elt0  +  y0OUifi  Sec  (*1  ■  *2^ 


1/2 


We  obtain  a  and  Bn  by  solving  the  last  equation  simultaneously  with 


60  "  a0  =  Vl^ 


2  2  2  2 

(which  is  equivalent  to  3^  -  +  | 6  -  | ) ,  once  for  a0  =  20  “  3q  and 
once  for  3Q  =  Cg  -  oig  (I  shall  not  write  down  the  derivation).  This 
quadrature  is  redundant,  practically  speaking,  since  the  most  feasible 
solution  seems  to  be  solving  A,  B,  a,  and  3,  as  stated  previously,  and 
then  obtaining  tp ^ ^ »  8*  0^ ,  and  Bg,  successively. 


If  we  now  consider  a  different  problem  and  prescribe  independently 
of  the  e  medium 


a  =  0  , 


3  •>  8  sin  9q  , 

fi  ■*■  1  , 

for  an  infinite  solution  plane  (where  0g  is  the  real  incidence  angle 
from  vertical),  then  the  usual  Fresnel  equations  determine  the  transmit¬ 
ted  amplitudes.  We  can  solve  the  0^  equations  given  above  by  way  of 

"iYo 

c°s  01  =  pQe 

(using  also  sin  0^)  for  Pg  and  y g.  Thus,  in  this  case  we  obtain 
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A  = 


P0U0 


1/2 


u(  n 


COS  yq  - 


<  sin  y„) 


B 


1/2, 
'0P0 


P„U„'  J'w(ic  cos  y  +  n  sin  y„) 


(59) 


and  this  solution  (trivially  allowing  for  a  different  e  la)t  convention) 

is  identical  to  that  given  by  Stratton,1®  for  the  well-known  problem 

just  stated.  (My  B,  A,  n,  x,  y  ,  and  8  correspond  to  Stratton's 
1/2  1/2 

p,  q,  (x^/uVq  ,  B.jAni  ,  Y/  and  a2>)  This  problem  is  a  homo¬ 
geneous  plane  wave  (ft  =  1)  for  which  the  limit  =  0  yields  Yg  =  *  = 
B  =  0,  A  *  0,  and  real  refraction  angle  i|i  =  0^.  In  this  case,  reac¬ 
tive  damping  (A=0,  B  *  0,  +  0)  is  not  possible. 


Returning  now  to  the  transmission-line  surface-wave  problem,  equa¬ 
tions  (59)  for  A  and  B  remain  true  but  a  and  6  are  no  longer  pre¬ 
scribed,  and  ft  t  1  has  been  restored.  This  exact  special  case  in  which 
e2  =  0  in  the  refraction  formalism  requires  that  A  now  vanish;  hence, 
we  require  that 


n  cos  Yg  ~  K  sin  Yg  =  0 


tan  Yg  =  n/K  *  °° 


1  .  S.  Stratton,  Electromagnetic  Theory,  McGraw  Hill,  Inc.  (1941), 

p  502  . 
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as  <  +  0;  hence,  a  =  0,  =  n/2,  A  =  0,  and 


1/2 .  ,  1  1/2 


p0  =  B0/U0  “(£1Q) 


cos  ei  =  -iB0/w(p0E1fi)  , 

sin  81  =  0/w(yoE , 


tan  0,  =  -8/ iB  =  E  /E 
1  0  -x  -z 


where  Bg  is  B  =  Bg  given  by  (33).  In  this  case  iJi  =  n/2 
and  =  0.  Defining  the  complex  refraction  angle  as 


61  =  -i9r  +  2  ' 


where  0^  is  a  real  angle,  one  finds  sin  01  =  cosh 
i  sinh  0r,  and 


0f  =  -coth  1  (8/Bq)  ; 


2  2  2  2  2 
hence,  cos  01  +  sin  0^  =  1  is  just  8  ”  Bq  =  poeiW 


18' 


because  1  =  | 1 1 . 


(pure  forward) 


0  ,  cos  0,  = 

r  1 


true  under 
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I  shall  digress  for  just  a  moment  to  give  a  general  proof  of  the 
absolute  condition  in  the  dispersion-relation  solution.  The  claim  is 
that,  in  general  (for  any  e^),  equation  (41)  is 


Re(cos^  01  +  sin^  9^)  =  1 


(60) 


which  is  identical  to 


2  2 

|  cos  9  +  sin  9^1  =  1  • 


(61) 


The 


proof  is  immediate  from  (46)  and  (47)  with  cos  0^ 


-n  +  in  , 
rx  lx 


sin  91  =  nrz  -  inIz.  We  can  verify  easily  that  (60)  is  identical  to 

2 

(61)  because  fn  nT  +  n  n„  )  =  0  as  a  result  of  using  (43)  and  (44) 

v  rz  Iz  rx  lx' 

in  elimination.  Thus,  for  ei  >0,  the  absolute  elimination  of  (43) 
means  only  that  ft  is  positive. 


To  obtain  SI  and  finish  the  lossless  solution  we  must  avoid  taking 


the  limit  +  0  directly  in  S2  (57).  The  convenient  additional  equation 
-*•  2 

is  the  |k|  norm  from  equation  (58)  or 


„2  2  2-1 

6  +  B0  =  a 


(62) 


2  2 

Thus,  since  |8  -  BQ| 


PgE.jU)  U,  one  has 
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=  V,f3 


B 


2.1/2 

O1 


On  expressing  the  left-hand  side  of  the  last  equation  exactly  as  was 
first  given  above  (between  eq  (59)  and  (60))  and  then  eliminating  Bq 
using  (62),  we  find 


1 26"  -  '  I"  "  = 


2  2-1  2 

Solving  uniquely  -26  +  ft  =  p^e^io  SI  one  obtains 


ft  =  -62/pQe1 0)2  +  [l  +  e4/(u0 


2  -.2  -i  1/2 


(63) 


(since  ft  is  real  and  the  plus  sign  is  unique  to  ft  positive  definite). 

1  /2 

In  the  |X^|  negligible  high-frequency  regime,  6  -  w(uoeo)  and  equa¬ 
tion  (63)  becomes 


ft  “  1  -  (e0/e^  +  (eo/eJ2/2  •  (64) 

The  ft  of  equation  (64)  is  close  to  1  for  real  soil  having  eC)/e1  -  1/6  or 
so  at  high  frequency. 

In  real  soil,  however,  high-frequency  a  and  6  can  have  a  slight 
effect  on  the  soil  fields  (as  if  a  -  0)  yet  a  can  be  still  appreciable 
unless  C2//ci  is  sufficiently  small.  In  order  to  see  what  constitutes 
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small  £ 2/z  t^le  following  high-frequency  approximations  to  Rg  and  Xg 
are  useful : 


,1/2 


>  1 


e2/£1 


<  1, 


ei/£o 


>  5,  (65) 


X 


high 


(66) 


The  gj/c  >  1  stipulation  is  essential  so  that  B  is  sufficiently  large. 
These  rough  estimates  are  good  to  about  20  percent  or  better  if 


(3/2)(e1/e0)"3/2w_1(u0/e0)1/2(^e}"1  <  0.1  .  (67) 

At  high  frequencies,  these  estimates  improve  to  a  few  percentage  points 

(factor  (68)  should  be  used  in  Xg).  As  frequency  increases,  C2^t‘\  be“ 

comes  smaller.11  Hence,  g  =  w/c,  and  a  =  (w/c)  sin  fR  /2w£  1.  The 

^  g  e ' 

estimates  come  from  equation  (44),  dropping  a  and  g,  and  using  A2  +  B2 
2 

=  u  q£  ^  Co  (where  Q  =  1).  This  yields 


\igh  =  °’5  a^0/E1^ 


1/2 


Bhigh  =  ^Vl)172'  lEl  “  £1 


^C.  L .  Longmire  and  K,  S,  Smith,  A  Universal  Impedance  for  Soils, 

Defense  Nuclear  Agency,  Topical  Report  DNA-3788T  (October  1975), 

1 

i 
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thus  equations  (65)  and  (66)  can  be  obtained  from  equations  (30)  and 


(31).  Incidentally,  taking  =  0,  tan  ij/^  =  6/A,  8  -  oj/c,  and  using 

^  for  A,  we  find  directly  from  (57)  that  ft  =  ( 1  +  4  Eg/e^) 

-1/2 

=  1  +  e„/2e  ,  in  agreement  with  ft  from  (64)  to  leading  order 

in  Eg/e^.  Since  (66)  is  entirely  independent  of  e^,  Xg  should  come 

from  -wu£,X  =  In  (33)  we  should  retain  the  l  term,  however, 

1  g  0  e 

2 

as  6  is  not  quite  negligible  with  respect  to  calculating  Xg.  The 
result  is  simply  an  ad. itional  factor  of 

h  -  v/£i)1/2  * 1  “  e0/2ei  (68) 

on  the  right-hand  side  of  equation  (66).  Finally,  with  (67)  well  satis¬ 
fied,  Rg  is  essentially  entirely  due  to  ohmic  loss.  The  tendency  at 
high  frequency  as  e2^ei  di-minislles  is  for  01  to  be  Purely  ohmic  while  B 
is  purely  reactive.  The  surprising  feature  of  real  soil  is  that  o  can 
be  appreciable  and  yet  e2/ei  can  be  small  in  a  "lossless"  situation. 

As  for  low-frequency  approximations,  we  again  drop  a  and  8  in  (44), 
taking  A  =  B.  The  familiar  result  for  Rg  =  Xg  is 

R.  =  X  =  (y  u)/2o)1/2w_1  .  (69) 

low  low  ^  0  ' 

The  approximation  can  be  good  (better  than  ~10  percent)  only  for  e2//£i  j 

10. 


- - ^  T — 
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4.  PURELY  VECTOR-POTENTIAL  MODEL  OF  TM  TRANSMISSION  LINE 


It  is  of  considerable  qualitative  interest  to  obtain  the  section  2 
and  soil  field  model  from  a  purely  vector-potential  TM  formulation  of 
the  wire-over-ground  problem.  In  our  region  of  macroscopic  electro- 
dynamics,  V»D  =  0  everywhere;  hence,  no  scalar  potential  need  be  intro¬ 
duced  explicitly  (it  is  zero). 


Because  V*A  =  0,  we  first  consider  a  TM  vector  potential  A  in 

region  1  (air)  volume  that  is  purely  radial.  The  radial  electric  field 
►  « 

E*r  =  E  (where  rr  is  outward  from  the  wire)  is 
-r 


-iwA 


rl 


cv 

—I —  =  E 
2ire„r  -r 


(70) 


where 


-ikzz  ^ 

the  wave  dependence  is  e  e  ,  implicit  in  V  of  equation  (15a), 


[  (x  -  h)2  +  y2] 1//2  ,  and 


A  is  the  incident  vector  potential  due  solely  to  the  current  I  on 


the  wire. 
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The  boundary  condition  at  the  wire  surface  is  satisfied  by  equation 
(70),  which  yields  normal  D  or  free  CV  waves  at  the  wire;  V  is 
specified  as  in  section  2.  This  must  be  so  since  =  0  and 
(curl  H)^  =  0  in  the  wire  (TM);  hence,  E  at  the  wire  surface  jumps 
abruptly  to  an  exceedingly  small  value  inside  the  metal. 

The  reaction  A  at  x  =  0  must  be  added  to  equation  (70).  Because  the 
reaction  must  instantaneously  yield  only  at  x  =  0  and  every  y  point, 
the  reflection  or  reaction  symmetry  of  the  divergenceless  axial  vector 

A 

A^r  must  be  opposite  that  of  an  optical  mirror  reflection  of  the  vector 
arrow,  so  that  y  projections  cancel  in  the  sum  and  Ey  =  0  identically 
over  x  =  0.  Then  the  normal  projection  of  equation  (70)  is  doubled 
at  x  =  0 : 


A 

-x 


1 


(71) 


2  2  1  /2 

where  rQ  =  (h  +  y  )  .  If  we  integrate  normal  over  y  at  x  =  0, 

from  equation  (71),  we  obtain 


ra 

-CV  =  /_ 


2he„iu)A  ,  dy 
0  -r  1 


which  is  exactly  satisfied  on  inserting  (70). 


(72) 
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The  integrand  of  (72)  is  a  transverse  magnetostatic  y  distribution 


so  that 


2x(0,y) 


2h(-CV) 
2n(h2  +  y]2 


(73) 


Hence,  H^(0,y)  has  a  distribution  like  that  of  (73),  with  -CV  replaced 
by  -I.  To  obtain  this  result,  we  substitute  -CV  =  -k  I/w  into  equation 
(73),  and  use  D^(0,y)  =  kzH^(0,y)/w  from  (12).  We  can  use  equation 
(73)  to  verify  (10),  since  0.5D  (0,0)w  =  -CV,  where  0.5D  (0,0)  is  the 

X  —  —X 

—  1  ,oo 

"conserved  mean  planar  field,”  that  is,  w  J  ^dy  Dx(0,y).  A  similar 

mean  value  definition  applies  to  and  also.  The  consistency 

check  here  is  that  in  the  contact  limit  h  =  a,  the  following  conditions 

hold:  w  =  2ira  and  D^(0,0)  =  -2CV/w,  as  shown  by  equation  (73);  thus, 

Dx(0,0)/2  =  -CV/w  is  the  negative  (r*x  =  -1)  of  the  radial  w  uniform 

D  around  the  wire  (D  is  also  the  mean  value).  The  line  wave  connec- 
-r  -r 

tion  to  the  region  2  soil  fields  of  section  2  is  now  complete  from  the 
air-side  approach  to  x  =  0.  One  can  obtain  all  field  amplitudes  from 


H  (described  in  (20)  and  below  (11)), 

-yg 


since  we  find  from  equation 


(10),  equation  (73),  and  H  (0,0)  =  wD  (0,0)/k 

-y  -x  z 


H  =  we  k_1 [0.5E  .(0,0)1  =  -wk_1CVw' 
-yg  0  z  L  — x 1  J  z  - 


where  CgE^^O.O)  =  (0,0) .  The  essential  explicit  minus  sign  of 

results  in  a  Poynting  flux  density  ( 0 ,y )H* ( 0 ,y )  into  the  ground. 
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Let  us  now  show  that  l  ,  Z .  ,  and  Z  are  already  included  in  the  TM 

e  l  g 

model  in  a  "consistently  exact"  way.  At  the  wire  and  conductor  sur¬ 
faces,  two  kinds  of  z-component  vector-potential  surface  A^  are  pres¬ 
ent:  (1)  A  „ ,  which  is  the  l  I  contribution,  and  (2)  A  for  which 

-zs  1  e-  -zs  2 

-iuiA  is  the  surface  electric  field  in  the  surface-impedance  boundary 

“ZS2 

condition . 


At  the  wire  surface  we  have  surface  H  =  where 


3  A  3  A 
-r 1  -zs2 


0-i  3  z 


3r 


r  >  a 


(74a) 


and 


Vi  = 


3  A  , 
-zs  1 

3r 


(74b) 


We  may  obtain  our  results  from  equations  (74a)  and  (74b)  with  r  locally 
free  and  then  pass  to  r  =  a.  The  radial  A^  derivative  from  equation 
(70)  is 


—iw 


3A 
-r  1 

3z 


2"e0r 


f-i0)i  I  -  Z.  I  j 
11  e-  i-' 


r  +  a 


(75) 


where  the  parentheses  in  equation  (75)  contain  the  part  (linear  super¬ 
position)  of  3V/3z  that  contributes  locally  at  the  wire  surface.  The 
first  term  in  equation  (75),  along  with  (70),  contributes  where 
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H  . 
-l 


I 

2ur 


r  =  a  . 


As  a  Taylor  expansion  shows,  A^  cannot  change  H  locally--nor  hence 
anywhere  else — inside  an  air-medium  pillbox  at  the  surface  that 
collapses  to  zero  radial  thickness.  The  second  term  in  equation  (75) 
thus  must  cancel  -iu)(  3 Azs2/^r )  from  equation  (74a).  This  yields 


E  =  -iwA 

-zwire  -zs2 


Z.l(-2nr  +  D) 

l- 


f 

g 


r  =  a  , 


(76) 


where  C  =  2ire^/f^.  The  correct  result  (independent  of  f^)  in  equation 

(76)  follows  if  constant  D  is  adjusted  so  that  f  =  £n(2h/a).  The 

g 

remaining  flzs1  must  necessarily  satisfy  (74b).  This  result,  where 
is  as  given  above  and  constant  D  is  as  in  (76),  gives 


A 

-zs 


1 


l  I 
e- 


(77) 


where  l  =  \i  f  /2tt.  Clearly,  A  has  an  ignorable  (zero)  local  induc- 
e  0  g  -zsl 

tion  contribution  to  Z  in  the  pillbox  argument,  and  the  sum  of  equa¬ 
tions  (74a)  and  (74b)  is  obviously  not  meaningful  (H  .  *  2H.). 

-wire  -l 


At  the  ground  surface  x  =  0,  the  surface  vector  potential  must 

be  A  =  0;  this  condition  is  required  so  that  local  A  yields  a 
— zs i  —zsl 

zero  contribution  to  equation  (77)  from  a  surface  integral  over  x  =  0. 
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Similarly,  at  x  =  0  we  have 


IJ 


„H 

0-y 


3A  „ 
-2s2 

3x 


3A 

-x 

3z 


(78a) 


and 


n0V°'Y>  - 


3A 
-zs  1 

3x 


(78b) 


The  A^  derivative  is,  in  general,  for  arbitrary  y, 


3  A 


-XU)  — —  =  - —  I  -  Z  i)  . 

3z  2  e-  g-' 

ne„r 
0  0 


(79) 


Equation  (79)  follows  from  A^  above,  (71),  (70)  and  3v/3z  acting  local¬ 
ly  at  x  *  0. 


The  3a^/3z  contribution  from  the  second  right-hand  term  of  equation 


(79)  must  cancel  with  -f3A  „/3x]  in  (78a).  Actually,  our  inclusion  of 

v  -zs2 

the  cancelling  terms  in  the  equations  is  not  necessary,  since  only  the 
one-line  boundary  condition  (Az  in  eq  (23))  at  x  =  0,  y  =  0  must 


result.  The  "cancellation  at  a  point"  condition  results  in  ‘-ZS2Y^X  = 

-f  2irf  lH  (0,y)Z  from  the  air  side,  where  E  _  =  -iwA  „  and  y  0  are 
v  g;  — y  rr  -zs2  -zs2 


implied.  Assuming  that  H  (0,y)z  =  E  (0,y)w  and  solving  the  last 

-y  g  -z 

derivative  equation  radially  by  way  of  3/3x  =  (-h/r)(3/3r)  on  SZS2^r^' 


where  E  (0,y)  =  E  (0,0)h2/r2, 

“Z  “Z 


we  find 


U nr  +  D) 


(80) 


E  (0,0) 

E  (r)  =  — - 

-zs2  f 

g 

Thus  equation  (80)  is  the  identity  E  (0,0)  as  r  r(0,0)  =  h,  where 

-z 

D  =  lu2  -  Jlna.  The  remaining  !) A^/3z  contribution  from  equation  (79) 

yields  g^H^(0,y),  thus  entirely  satifying  equation  (78a). 


Analogous  to  (74b),  equation  (78b)  simply  defines  the  local  normal 

derivative  of  A  .  A  is  again  ignorable. 

"Zs I  -zs  i 


It  is  obvious  that  f  =  Jtn(2h/a)  is  not  required  for  model  consist- 

g 

ency.  In  fact  we  must  replace  this  l  with  the  resulting  from  the 
more  exact 


The  origin  of  air  volume  Z^  can  be  attributed  to  a  linearly  super¬ 
posed  vector  potential  A^z  that  contributes  the  surface  Azgi  consid¬ 
ered  previously.  Let  subscript  t  denote  transverse  plane  vector 

2  2  2  2 

components.  Then  V  =  +  [  3  /3z  ) .  The  vector  potential  A  ^ 

-ikzz 

carries  the  suppressed  e  e  dependence.  A  satisfies  the 

—zi 

everywhere  instantaneous  Faraday  Law 

vtxt"iMS*i*)  =  ~iu>vo?ot  ' 
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where  H  =  [curl(  A^z)  ]fc.  Evidently  the  wave  equation  (V  +  k^F  = 

0  =  V  F  is  satisfied  for  F  =  Az-)<HQt  and  w^ere  Eot  = 

-iwA  _ z .  The  last  example  shows  that  A  .z  is  a  source-free, 

-z 1  -z  1 

ignorable  divergence  axial  vector;  hence,  the  "incident  plus  reaction" 
boundary  condition  Azsi  =  ®  necessarily  is  required  at  x  =  0  in  solv¬ 
ing  V^A  „  =  0.  We  may  note  that  A  _  creates  zero  H  ;  hence  V«H  = 
t-z 1  — z 1  -z  -0 

-*■  -*• 

V  *H  =  0.  Thus,  V*H  =  0  is  satisfied  completely  in  the  TM  model. 
The  meaning  of  =  0  is  that  3^A  ^/3x  3y  is  the  identity,  plac¬ 

ing  no  restriction  on  A  „ . 

-z  1 


A  useful  boundary  relation  at  wire  surface  is 


leI  a[curl(A2lz)]*t 
__  =  - 

9  9 


(82) 


where  t  is  the  tangent  vector  along  a  wire  cross-section  perimeter  (A^ 
is  TM  uniform  in  <j>  at  the  wire  surface).  A  numerical  solution  is  not 


necessary,  because  the  standard  circular  equipotential  approximate 

2 

solution  (vector-potential  version1'1)  of  v  A^  =  0  here  results  in  a 

surface  potential  value  A  .  =  D  f  constant  on  the  circle  of  wire 

— z 1  -0  g 

surface  ( f  is  defined  in  eq  (81)).  We  may  construct  a  locally  La- 
placian  wire  coordinate-system  solution  for  r  =  a  of  A^  =  D^f^a/r, 
and  with  [curl(  A^z)  ]  ■*  -(SA^/ar),  the  boundary  relation  (82)  yields 
Dq  =  pQI/ 2it  as  r  ♦  a. 

l2See,  for  example,  S.  Ramo  and  J.  R,  Whinner y,  Fields  and  Waves  in 
Modern  Radio,  second  edition,  John  Wiley  and  Sons,  Inc.  (1953),  p  138. 
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alone 


Equations  (78b)  and  (74b),  where  H. .  .  =  H.,  reveal  that  H-. 

—  Ut.<J)  —1  —ot. 

-► 

must  yield  the  correct  H  in  the  air  volume  that  is  continuous  in  the 
tangential  component  at  the  wire  surface  and  at  the  x  =  0  surface,  on 
approach  to  the  surfaces  from  the  air.  Clearly,  then,  H^(r)<}>  = 
I4>/2TTr,  where  H^(r)  =  p^  ^  (  3A^_  ^/Sz )  .  This  result  is  nonzero  only  at 
r  =  a  and  r  =  r^.  The  total  H,  at  x  =  0  is  2H^(rp),  and  H^(0,y)  is  in¬ 
deed  the  y  projection  2tr(rg)<t>*y  =  p^  (SA^/dz)  ^ ,  where  <j>»y  =  -h/r^, 
and  where  (SA^/Sz)^  comes  entirely  from  the  noncancelled  first  right- 

hand  term  of  (79).  That  H.  and  hence  A  „  must  vanish  in  the  air 

-l  -r  1 

volume  follows  from  curl  curlfH.ifcl  =  32H.d>/3z2  =  U)2p.e„H.d>  and  curl 

-it  0  O-i 

curlfA^r)  =  -32Ar1r/3z2  =  0.  The  propagating  kz  dependence  of  is 
essential . 

In  the  TM  vector-potential  model,  A^,  as  well  as  3Ar1/3z,  is  actu¬ 
ally  zero  in  the  air  volume,  with  A^  in  normal  D  also  entering  the 
model  only  as  wire  and  x  =  0  surface  values.  The  air  volume  vanishing 
of  Ar1  is  absolutely  essential  physically;  otherwise,  at  x  =  0  +  e  (e 

is  infinitesimal)  there  is  an  additional  (to  eq  (73))  D  contribution, 

-x 

2  2 

k  f -i<oA  .  ]/p„w  where  A  „  =  -hA  ,/r.  in  3A  „/3z.  This  false  D  ,  which 
zK  -x1;  0  — xl  -r 1  0  — x 1  -x 

is  made  continuous  (suppose  superposition)  with  a  from  the  soil, 

requires  common  air-soil  k^  =  uifp^Eg)  .  The  surface  curl, 

fcurl  1  =  io)D  (0,y),  is  continuous  in  3H  /3z  at  x  =  0  when 

v  -Ot'x  -x  -y 

H„  =  H  (0,y).  The  remaining  x  =  0  curls  of  are  fcurl  H_  1"  =  0 
-Ot  -y  ’  -Ot  -0t;z 

(where  A  .  =  0)  and  fcurl  H„  )  =  0  (normal  H„  vanishes).  At  the 

-zsl  v  -Ot'y  -Ot 
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wire,  taking  H  ,  only,  we  similarly  find  fcurl  H„  1  continuous, 
-Ot*  1  -Ot'r 

(curl  H  )  =  0  and  (curl  H  )  =  ( 1/r ) ( 9/3r ) (rH  )  zero  with  local 

" Ut.  (p  “ U  u  Z  t  (I) 

Sot*  =  I/2llr '  r  =  a* 

In  the  vector-potential  model,  the  Neumann  analog  tangential  is 

redundantly  equal  to  the  Dirichlet  H  at  wire  surface  and  x  =  0.  We 

can  conclude  from  the  uniqueness  theorem  that  the  H  resulting  from 
2 

V  A  .  =  0  must  yield  the  identical  A  „  if  the  Dirichlet  is  pre- 
t-zl  -2 1  -Ot 

scribed  as  a  surface-tangential  H  condition,  thus  defining  the  Neumann 
problem.  The  Neumann  problem  then  becomes  a  straw  man  that  need  not  be 
considered. 

There  is,  however,  still  an  unresolved  question  concerning  the  exact 
Dx(0,y)  distribution  that  is  dependent  on  approximations  made  in  solv¬ 
ing  the  Dirichlet  A  ^  problem.  It  must  be  true  that  equation  (73)  is 
correct  since  the  traveling-wave  vector-potential  solution  is  physically 
unique  already.  A  complementary  (to  the  l ^  solution)  standard  approxi¬ 
mate  capacitance  solution12  can  be  devised  that  yields  what  we  call  the 
exact  capacitance  C,  based  on  a  circular  equipotential  at  the  wire 
surface;  however,  this  solution  cannot  predict  the  distribution  equation 
(73)  as  h/a  -*■  1,  since  both  static  formulations  result  in  equation  (73) 
with  h  replaced  by  d^  =  a[(h/a)^  -  and  w  =  2Tidg.  A  vanishing 

12See,  for  example,  S.  Ramo  and  J.  R.  Whinnery,  Fields  and  Waves  in 
Modern  Radio,  second  edition,  John  Wiley  and  Sons,  Inc.  (1953),  p  138. 
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w  as  h/a  ->-1  is,  of  course,  absurd.  One  has  to  give  up  h/a  -*■  1  to 
recover  d^  -*•  h  and  suppose  that  a  *  0  for  h/a  fixed  infinitesimal. 

Although  this  is  not  a  proof  of  the  exactness  of  equation  (73)  with 
respect  to  fq,  we  can  show  formally  in  a  special  case  that  of  (81) 
arises  from  equation  (73)  and  the  unit-voltaqe  line-charge  Green  func¬ 
tion  that  vanishes  at  the  observer's  location  where  the  source  is  at  the 
origin,  namely 


g 


where  p  is  the  observer  p,pQ  source  coordinate  with  respect  to  a  wire- 
based  coordinate  system.  We  calculate  the  (vanishing)  transverse  deri¬ 
vative  at  x  =  0 ,  y  =  0 , 

9G 

Dy(0,0)  =  -/^  dy'  Dx(0,y')  ,  (84) 

by  taking  equation  (73)  for  both  and  D^.  Taking  p  =  h,  (j>  =  0 ,  and 
3Gg/3y  =  VG^*y,  we  may  replace  the  last  equation  with  1  =  -/dp^ 
d<|>0p0  cos2  (f>Q  sin  <J>  aGQ/3 p t<5C4>0)P01  ]•  If  sin  is  substituted  for 
sin  $  and  coordinate  pQ  for  p  in  the  GQ  derivative,  the  result  is 
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fg  =  /; du^2 "  o 1/2 


0 


where  u  =  h/a;  hence,  f  is  the  limit  h/a  -*•  1  of  equation  (81),  because 

9 

of  £n( 1/1). 

The  TM  model  predicts  that  the  TE-related  amplitudes  of 

H^,  Ey,  and  are  exactly  zero.  A  formal  continuous  y-integrated 

matching  of  assumed  soil  and  at  x  =  0  (with  air  fields  arising 

from  Cartesian  components  due  solely  to  surface  A^)  results  in  soil 

field  amplitudes  H  and  H  being  dependent  solely  on  the  E^  amplitude 

(since  k  =0).  Thus,  the  formal  antisymmetry  in  y  of  H  and  H  in  the 
y  -x  -z 

air  from  A  ^  is  not  the  actual  reason  that  the  y-integrated  amplitudes 
vanish.  The  vanishing  is  due  to  the  stronger  condition  of  pointwise 
vanishing  required  of  E^.  The  vector-potential  A^  produces  zero 

because  of  the  vanishing  of  at  x  =  0. 

It  is  natural  that  the  TM  model  can  create  no  TE  waves,  since  the  y 
antisymmetry  in  E  and  H  is  in  fact  just  the  opposite  of  the  true 

y  x 

symmetry  in  E  and  H  at  x  =  0  that  would  arise  from  ground  reaction  of 

y  x 

the  circularly  uniform  A^  vector-potential  component.  This  symmetry 

yields  a  forward  Poynting  flux  down  the  line  from  E  and  H  .  Observe, 

y  x 

■f  in 

however,  that  if  we  examine  (curl  H),  inside  the  metal  wire,  k  =  0 

m  Z 

again  forces  Hf  to  jump  as  Ef  did  in  the  TM  case.  The  wire  surface 
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value  of 
hence, 


Hf  must  then  be  exceedingly  small.  H^r  is  an  axial  vector; 
from  such  an  H^_r  is  also  opposite  the  TM  symmetry  in  y  at  x  = 


0. 


At  the  same  time,  (curl  H)  in  the  metal  wire  cannot  yield  "zero"  E 

r  r 

unless  3H  /3 <p  is  "zero"  or  circularly  uniform  H  .  But  then  E  must, 

Z  Z  (J) 

from  (curl  H)  ,  also  have  the  circular  uniformity  of  H  at  th''  metal 

m  Z 

surface.  Now  E,  comes  from  vector  potential  A  ,  which  thus  contributes 

<P  <f> 

consistently  zero  divergence  to  V*A  =  0.  Consequently,  the  H  from  A 

z  9 

is  "zero"  because  <4  is  uniform. 


If  H  and  H  are  "zero,"  the  E,  is  also  "zero"  at  the  wire  surface; 

z  X  (f) 

hence,  (jj-uniform  TE  excitations  should  be  negligibly  produced  by  the 
wire.  Non-(j)-uniform  excitations  should  be  negligibly  excited  unier 
electrically  thin  wire  conditions,  in  the  absence  of  line  discontinui¬ 
ties  . 


It  is  peculiar  in  the  Wait  wire-over-ground  formulation2  that  the 
field  components  H^,  E^,  and  are  present,  yet  the  only  wire  boundary 

condition  is  that  E  =0  along  the  wire.  How  then  can  the  TM  excitation 

-z 

from  the  wire  excite  these  TE  components  by  way  of  ground  interaction? 
The  answer  to  the  question  is  that  the  Wait  formulation  should  not  allow 


^James  R.  Wait,  Radio  Science,  7,  No.  6  (June  1972),  675-679, 
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nonzero  H^,  E^,  and  H  .  We  can  find  no  error  in  the  following  proce¬ 
dure.  Take  the  y  derivatives  given  in  Wait's  paper2  under  the  sign  and 
integrate  fields  over  dy  from  -“  to  this  yields  a  6  function  to  be 
integrated  over  X,  that  is,  6(k^),  =  X  =  0.  Assume  that  P1  =  U2  = 

pg.  The  field  continuity  at  x  =  0  then  yields  the  result  M(0)  =  -N  ( 0 ) 

(from  three  equations  stating  the  H^,  H^,  and  E  continuity).  Hence,  it 

follows  immediately  that  M(0)  =  N(0)  =  0  and  that  H  ,  H  ,  and  E  are 

-x  -z  -y 

identically  zero  pointwise  over  x  =  0.  The  normal  D  continuity 
yields  (redundantly) 


2  2 

where  e,  =  e  .  and  e„  =  e  ...  and  0  is  the  negative  of  our  y  .  If  we 
-1  air  -2  soil  ^  ' 

3  3  2  2 

require  that  (e1  -  =  k(e1  -  £2)'  solution  of  the  last  equation 

2  2  2  2  2 
is  K  =  e,  =  e  and  0  =  e,  =  k„.  But  then  E  =  fk„  -  0  Itt  =  0 

-1-2  0-11  -z  v  1  ' - 

everywhere.  This  result  can  be  understood  in  the  implied  lossless  case 
(£1  =  e2>  since  a  lossless  TE  solution  has  zero  axial  electric  field. 

The  above  X  =0  modification  of  the  Wait  formulation  yields  (be¬ 
fore  we  set  e1  =  e^)  the  result  that 


2 James  R.  Wait,  Radio  Science,  7_,  No,  6  ( June  1972),  675-679. 
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0 


/I  dv  *z<°'y>  .(v o2£i  -  e2)1/2  [i  +  R<o>] 


Jl  dy  sy(°*yJ 


e  ^ u>[  1  -  R(  0 )  ] 


But  we  require  R(0)  =  0  (actually  1  -  R(0)  =  T(0)  +  1),  so  that 

,00 

H .  =  H„  at  x  =  0  results*  in  I  H,  dy  =  -I;  hence,  the  Hertz  poten- 

- 1y  -2y  -1y 

tial  ir  in  Wait's  formulation  can  result  formally  in  our  Z^,  except 

that  e  . ,  in  our  Z„  is  replaced  by  e  . 

soil  9  air 


A  final  remark  concerning  the  vector-potential  model:  the  additive 
ground  interaction  of  H  with  respect  to  the  incident  wire  H  will  not 
allow  what  is  sometimes  called  an  "antenna  mode"  with  co-directional 
wire  current  I  and  return  current  f  dy  dxfiweE  1.  So  long 

*  —oo  •  — oo  v  —  z ' 

as  k^  is  zero  in  the  air,  the  transmission-line  solution  seems  to  be  the 
solution  to  the  free-wave  wire-over-ground  problem. 


5.  CONCLUDING  REMARKS 


Although  there  is,  at  present,  a  lack  of  suitable  cw  data  for  a(oj) 
and  B(w)  concerning  wires  close  to  earth,  we  have  some  limited  evidence 
that  our  Zg  is  realistic.  Gray7  has  reported  propagation  velocities 
versus  wire  height  over  earth,  experimentally  deduced  from  timing  of 


7R.  F.  Gray,  Nuclear  Electromagnetic  Pulse  Simulation  by  Point  Source 
Injection  Techniques  for  Shielded  and  Unshielded  Penetrations ,  Harry 
Diamond  Laboratories,  HDL-TR-1737  (December  1975), 

*The  coefficient  of  the  Hertz  potentials  given  by  Wait  (ref  2)  fol¬ 
lowing  equation  (8)  is  incorrect  and  should  be  multipl ied  by  a  factor 
of  2  . 
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reflection  pulses  at  cable  center  on  an  open  40-ft-long  cable  of  radius 

1  cm.  Using  a  semi-empirical  £(u>),^  ignoring  dispersion,  and  taking 
7 

frequency  f  =  10  Hz,  we  find  an  oj/cfJ  of  0.59,  compared  to  measured  0.55 
at  h  =  0.05  m  with  a  soil  moisture  content  of  50  percent.  Our  oj/cB 
result  closely  follows  the  measured  result  upward  in  height  (<1  m)  to 
within  10  percent,  the  high-line  agreement  being  (always)  accurately  1. 


iJC.  L.  Longmire  and  K.  S.  Smith,  A  Universal  Impedance  for  Soils, 
Defense  Nuclear  Agency,  Topical  Report  DNA-3788T  (October  1975). 
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APPENDIX  A 


The  body  of  this  report  addresses  the  problem  of  calculating  the 
response  of  a  horizontal  wire  over  earth  to  incident  electromagnetic 
fields.  In  an  approximate  solution  of  this  problem,  the  wire  over  earth 
is  regarded  as  a  transmission  line;  the  basic  unknown  in  this  solution 
is  the  relationship  between  the  line  waves  and  the  electromagnetic 
fields  in  the  earth:  that  is,  the  surface  impedance  of  the  ground 
plane. 

Equation  (3)  in  the  body  of  the  report  describes  Zg,  the  per-unit- 
length  surface  impedance  of  the  finite  conducting  earth.  In  this 
appendix  we  derive  equation  (3)  by  solving  equation  (21)  self- 
consistently  with  y2 »  where  y  is  the  line  propagation  constant  and  a 
function  of  Zg.  See  section  2  in  the  body  of  the  report  for  a  detailed 
description  of  the  model  and  definitions  of  symbols. 


Inserting  equation  (1)  into  equation  (21)  and  squaring  gives 


Z  +  c.e 

g  1 


Z  -  c  e 

g  2 


(A-1  ) 


The  parameters  c^  c^,  and  <p2  are  given  in  equations  (4)  through 

(10)  in  the  body  of  the  paper.  Basically,  we  can  solve  (A-1)  at 
the  point  =  ^2  =  ^  (this  is  a  physical  solution  for  and 

( l +  &e)c  =  UqE^)«  We  can  then  analytically  continue  the  resulting 
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appendix  a 


Z  to  the  desired  parameters, 

g 

of  Zg  is  physically  correct. 


having  verified  that  the  e2/,£1  >> 


limit 


Proceeding  algebraically,  with 
(A-1)  and  find 


<*>  =  o  understood,  we  solve 

¥2 


Z 

g 


i<t> 

0 . 5c  ^e 


(A-2) 


where 


^  =  4>2 


“  2*1 


( A-3) 


(The  positive  root  sign  is  correct,  as  later  shown.) 
gives 


Now  equation  ( A-2 ) 


Z  =  -0.5c,  cos 

g  1 


+  0.50^  cos  [<t>1  +  0,5  tan 


(r)  ] 


( A-4) 


+  i(-0.5c1  sin 


+  0.50^  tin 


[♦, 


+0.5  tan 


(r)]) 


with 


[(1 


+  4c2  cos 


2  .  2 
+  16c2  sm 


♦/«;] 


4 1 1/4 
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and 


r 


4c„  sin 
2 


w{ 


c1  +  4c2  cos 


♦  ) 


From  (A-4)  we  verify  that 


■  */2)  =  XgCt,)  ' 


which  is  the  identical  phase-following  behavior  of 


Zg  =  lzgl|c°s  [  4> ,  +  +  1  sin  [♦1  +  5(<t>2)]|  •  (A-5) 

At  <(>1  =  0  we  can  determine  the  functional  form  6  ( )  of  the  phase  shift 
6 ,  if  we  require  that 


5  + 


0, 


X  =  0 

g 


(A-6) 


From  (A-6)  and  ( A-3 )  we  find  immediately 


6  =  -0.5<J>2 


(  A- 7 ) 


Furthermore,  with  cos  ($)  =  1  (where  =  0)  we  have  from  (A-2),  just 


|Zg|  =  0 • 5c1  [  - 1  +  (1  +  4c2/c^)1/2] 


( A-8 ) 
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which,  together  with  equations  (A-5)  and  (A-7),  yields  the  general 
solution  (3)  given  in  the  text. 

2 

If  e2/e1  >>  1/  4>1  =  <j>2  =  n/2,  and  c2  =  wi^/w  oQ  (where  c2/c1  >>  1), 
we  see  that  equation  (3)  gives  (with  negligible  error)  just 

wZ  -*•  (nfy  /o  )1,/2(  1  +  i)  ; 
g  v  0  0 

this  is  the  familiar  textbook  one-dimensional  result  which  can  also  be 

2 

obtained  (to  the  same  approximation)  by  simply  dropping  y  in  equation 


APPENDIX  B 


In  the  body  of  the  text  we  showed  that  the  electromagnetic  fields  in 

the  earth  er  inhomogeneous  plane  waves  with  wave  vector  k  = 

+  +2  ■>2  2  _ 

-  iaQn.  The  square  k  factors  to  k  =  where  is  the 

+2 

unsigned  ft  such  that  ft  =  |ft  |  .  The  imaginary  part  of  the  k  dispersion 
♦  ■¥  2  (  \ 

relation  is  -2a^30m»ni  =  this  requires  that  ft  -►  ft  since, 

2  2 

for  example,  3Q  <  aQ  and  ft^  <  0 ,  because  A  <  B  in  the  metal  limit 

-*•2  2 

where  a  »  3  =  0.  We  showed  that  k  =  Pgio  eft  has  the  correct  solution 
-*•2  -*• 2  2-*- 2 

because  n  =  1  in  k  =  k  n  ;  this  results  in  the  other  half  of  the 

2  2  2 

dispersion  relation  solution  being  3g  -  a  =  e^ft.  If  we  take  the 

2  2 

absolute  value  of  3.  -  a„,  this  last  equality  allows  ft  >  ft  to  have 

0  0  u 

+2  2„ 

the  same  consequence  as  the  real  part  of  k  =  pQu)  ft^  taken  as  an  abso¬ 
lute-value  equality. 


In  order  to  solve  explicitly  for  ft,  however,  it  is  necessary  to 
observe  that  equation  (52)  holds;  also,  equation  (50)  follows  from  the 
Poynting  divergence  relation. 


To  derive  equation  (52)  we  note  that  the  factorization  of  ft, 


u 

-1 


♦2 

in  k  changes  no  physical  result.  Thus  we  must  first  factor  out  ft 

'u 

■*2 

from  k  in  the  equation  of  squared  amplitudes  that  is  formed  from  the 


squared  Maxwell  curl  E  equation,  namely 


-1+2+2 
a  k  E 

S2  =  _y _ z_ 

2  2 

Mo“ 


(B-1  ) 


67 


APPENDIX  B 


so  that  ft  Vft  eu„ in  equation  ( B—  1 )  gives 
u  ^  u  0  1 


*2 

H  = 


tl2 


( B-2 ) 


Equation  (B-2)  also  results  from  the  squared  Maxwell  curl  H  equation 

with  factored  Si  ^k2  aqain  substituted  for  k2.  In  the  homogeneous  wave 
u 

+  +  >2  2  . 
equation  solutions  for  E  and  H,  the  vanishing  factor  k  -  PqEid  becomes 

-1*2  2  *2  2 
SI  k  u„ea)  =  0  if  k  =  |j„eoj  ft  .  It  is  then  true  that 
u  0  0  u 


I  E2  |  =  ft  I E  | 2 


( B-3 ) 


since  substituting  equation  (B-3)  into  the  absolute  value  of  either 
equation  (B-1)  or  equation  (B-2)  yields  the  same  factorization- 
independent  result,  namely, 


ifl2  («w0)‘ 


I  H  |  2 


Ik2! 


which  is  the  first  equality  in  equation  (52);  |H 


>2, 


I  H  |  2  . 


From  the  two  Maxwell  equations  and  equation  (B-3)  we  also  find 


ft"1(u2u)2|H2|/ft"1|k2|)  =  |k  | 2  |h  |  2/oj2  |  e  | 


1 
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in  which  the  curl  H  equation  (curl  E  equation)  is  used  in  the  right 
(left)  side  of  the  last  equation.  Thus,  the  consequence  is 


I*2  I 


Ik|2 


.2  2 
1 


( B-4 ) 


which  is  truly  independent  of  ft  from  the  vector  norms  of  equation  (58). 
We  then  use  the  n  ,n  formalism  to  evaluate  (B-4),  i.e.. 


*  + 
k  =  k 


-  ikI  , 


t  1/2  r  ■*  \ 

kr  =  uQ  o)[nnr  -  <n j J  , 


i  1/2  (  +  +  \ 

kj  =  W0  +  nnrJ 


where  n^  and  n^  are  equations  (46)  and  (47),  respectively.  We  then  find 


>2  ^  >2  -2 
n  +  n,  =  ft 
r  I 


( B-5 ) 


Further  evaluating  equation  (B-5),  we  find 


1  =  (a2  +  B2  +  B2  +  a2)(n2  +  <2)/p0ai2|e|2  . 


( B-6 ) 


This  result  proves,  in  fact,  our  assertion  (52),  since  the  norm  is 
unique  to  equations  (B-6)  and  (58).  We  saw  from  the  beginning,  in  the 
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complex  index  of  refraction  development,  that  ft  cancels  out  formally, 

-► 

giving  an  algebraic  identity  between  k  forms. 

The  point  is  that  equation  (B-5)  must  evidently  be  true 

-2 

independently;  therefore  ft  factors  out  in  (B-5)  leaving  the  remaining 

2  2 

explicit  ft  dependence  entirely  in  n  +  k  =  |e|ft  in  equation  (B-7). 
-► 

The  squared  E  norm  relation  above  shows  that  ft  =  1  can  only  be  true 
if  E  has  one  vector  component. 


70 


DISTRIBUTION 


ADMINISTRATOR 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTN  DTIC-DDA  (12  COPIES) 

CAMERON  STATION,  BUILDING  5 
ALEXANDRIA,  VA  22314 

COMMANDER 

US  ARMY  RSCH  &  STD  GP  (EUR) 

ATTN  CHIEF,  PHYSICS  &  MATH  BRANCH 
FPO  NEW  YORK  09510 

COMMANDER 

US  ARMY  MISSILE  &  MUNITIONS 
CENTER  &  SCHOOL 
ATTN  ATSK-CTD-F 
REDSTONE  ARSENAL,  AL  35809 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS 
ACTIVITY 
ATTN  DRXSY-MP 
ATTN  DRXSY-PO 

ABERDEEN  PROVING  GROUND,  MD  21005 
DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH 
LABORATORY 

ATTN  DRDAR-TSB-S  (STINFO) 

ATTN  DRXBR-AM,  W.  VANANTWERP 
ATTN  DRSTE-EL 
ATTN  DRDAR-BLE 

ABERDEEN  PROVING  GROUND,  MD  21005 

US  ARMY  ELECTRONICS  TECHNOLOGY 
&  DEVICES  LA  ORATORY 
ATTN  DELET-DD 
ATTN  DELSD-L 
FT  MONMOUTH,  NJ  07703 

HQ,  USAF/SAMI 
WASHINGTON,  DC  20330 

TELEDYNE  BROWN  ENGINEERING 
CUMMINGS  RESEARCH  PARK 
ATTN  DR.  MELVIN  L.  PRICE,  MS-44 
HUNTSVILLE,  AL  35807 

BOEING  MILITARY  AIRPLANE  CO 
ATTN  A.  LAWRENCE  WELLER 
M/S:  K75-50 

3801  SOUTH  OLIVER 
WICHITA,  KS  67210 

ENGINEERING  SOCIETIES  LIBRARY 
345  EAST  47TH  STREET 
ATTN  ACQUISITIONS  DEPARTMENT 
NEW  YORK,  NY  10017 


DIRECTOR 

ARMED  FORCES  RADIOBIOLOGY 
RESEARCH  INSTITUTE 
DEFENSE  NUCLEAR  AGENCY 
NATIONAL  NAVAL  MEDICAL  CENTER 
ATTN  RESEARCH  PROGRAM 

COORDINATING  OFFICER 
BETHESDA,  MD  20014 

ASSISTANT  TO  THE  SECRETARY  OF  DEFENSE 
ATOMIC  ENERGY 
ATTN  EXECUTIVE  ASSISTANT 
WASHINGTON,  DC  20301 

DIRECTOR 

DEFENSE  ADVANCED  RSCH  PROJ  AGENCY 
ARCHITECT  BUILDING 
1400  WILSON  BLVD 
ATTN  TIO 

ARLINGTON,  VA  22209 

FEDERAL  EMERGENCY  MANAGEMENT  AGENCY 
ATTN  JAMES  W.  KERR 

MITIGATION  &  RESEARCH 
WASHINGTON,  DC  20472 

DIRECTOR 

DEFENSE  COMMUNICATIONS  AGENCY 
ATTN  CODE  C312 
ATTN  CODE  C313 
WASHINGTON,  DC  20305 

DEFENSE  COMMUNICATIONS  ENGINEERING 
CENTER 

1860  WIEHLE  AVENUE 

ATTN  CODE  R720,  C.  STANSBERRY 

ATTN  CODE  R123,  RSCH  LIB 

ATTN  CODE  R400 

RESTON,  VA  22090 

DIRECTOR 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN  RDS-3A 

ATTN  RDS-3A4,  POMPONIO  PLAZA 
WASHINGTON,  DC  20301 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
ATTN  DDST,  E.  E.  CONRAD,  DEP  DIR 
SCI  &  TECH 

ATTN  RAEV,  ELECTRONICS  VULNERABILITY 
DIV 

ATTN  TITL,  TECHNICAL  LIBRARY  DIV 
ATTN  RAEE,  EMP  EFFECTS  DIV 
WASHINGTON,  DC  20305 


71 


DISTRIBUTION  (Cont'd) 


COMMANDER 
FIELD  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
ATTN  FCPR 

ATTN  FCSPM,  J*  SMITH 
ATTN  FCLMC 

KIRKLAND  AFB,  NM  87115 
DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
ATTN  TTV 

KIRTLAND  AFB,  NM  87115 

JOINT  CHIEFS  OF  STAFF 
ATTN  J-3 

WASHINGTON,  DC  20301 
DIRECTOR 

JOINT  STRATEGIC  TARGET  PLANNING 
STAFF,  JCS 
ATTN  JSAS 
ATTN  JPST 

ATTN  NRI-STINFO  LIBRARY 
OFFUTT  AFB 
OMAHA,  NB  681 1 3 

NATIONAL  COMMUNICATIONS  SYSTEM 
OFFICE  OF  THE  MANAGER 
DEPARTMENT  OF  DEFENSE 
ATTN  NCS-TS,  CHARLES  D.  BODSON 
WASHINGTON,  DC  20305 

DIRECTOR 

NATIONAL  SECURITY  AGENCY 
DEPARTMENT  OF  DEFENSE 
ATTN  R-52,  0.  VAN  GUNTEN 
ATTN  S232,  D.  VINCENT 
FT  MEAD,  MD  20755 

UNDER  SECRETARY  OF  DEF  FOR  RSCH  &  ENGRG 
DEPARTMENT  OF  DEFENSE 
ATTN  G.  BARSE 
ATTN  S&SS  (OS) 

WASHINGTON,  DC  20301 

COMMANDER 
BMD  SYSTEM  COMMAND 
DEPARTMENT  OF  THE  ARMY 
PO  BOX  1500 
ATTN  BMDSC-AOLIB 
HUNTSVILLE,  AL  35807 

COMMANDER 

ERADCOM  TECHNICAL  SUPPORT  ACTIVITY 

DEPARTMENT  OF  THE  ARMY 

ATTN  DRDCO-COM-ME,  G.  GAULE 

ATTN  TECHNICAL  LIBRARY  DIV 

ATTN  DELCS-K,  A.  COHEN 

ATTN  DELET-IR,  E.  HUNTER 

FT  MONMOUTH,  NJ  07703 


COMMANDER 

US  ARMY  ARMOR  CENTER 
ATTN  TECHNICAL  LIBRARY 
FT  KNOX,  KY  40121 

COMMANDER 

US  ARMY  COMM-ELEC  ENGRG  INSTAL 
AGENCY 

ATTN  CCC-PRSO-S 

ATTN  CCC-CED-SES 

FT  HUACHUCA,  AZ  85613 

COMMANDER 

US  ARMY  COMMUNICATIONS  COMMAND 
COMBAT  DEVELOPMENT  DIVISION 
ATTN  ATS I-CD-MD 
FT  HUACHUCA,  AZ  8561 3 

CHIEF 

US  ARMY  COMMUNICATIONS  SYS  AGENCY 
ATTN  CCM-RD-T, 

ATTN  CCM-AD-SV 
FT  MONMOUTH,  NJ  07703 

PROJECT  OFFICER 
US  ARMY  COMMUNICATIONS  RES 
&  DEV  COMMAND 
ATTN  DRCPM-ATC 
ATTN  DRCPM-TDS-BSI 
FT  MONMOUTH,  NJ  07703 

DIVISION  ENGINEER 

US  ARMY  ENGINEER  DIV,  HUNTSVILLE 

PO  BOX  1600,  WEST  STATION 

ATTN  HNDED-SR 

ATTN  A.  T.  BOLT 

HUNTSVILLE,  AL  35807 

US  ARMY  INTEL  THREAT 
ANALYSIS  DETACHMENT 
ROOM  2201 ,  BLDG  A 
ARLINGTON  HALL  STATION 
ATTN  RM  2200,  BLDG  A 
ARLINGTON,  VA  22212 

COMMANDER 

US  ARMY  INTELLIGENCE  &  SEC  CMD 
ARLINGTON  HALL  STATION 
4000  ARLINGTON  BLVD 
ATTN  TECHNICAL  LIBRARY 
ATTN  TECH  INFO  FAC 
ARLINGTON,  VA  22212 

COMMANDER 

US  ARMY  MISSILE  COMMAND 

ATTN  DRCPM-PE-EA,  WALLACE  O.  WAGNER 

ATTN  DRCPM-PE-EG,  WILLIAM  B.  JOHNSON 

ATTN  DRDMI-TBD 

ATTN  DRDMI-EAA 

REDSTONE  ARSENAL,  AL  35809 


72 


DISTRIBUTION  (Cont’d) 


COMMANDER 

US  ARMY  TEST  S  EVALUATION  COMMAND 
ATTN  DRSTE-FA 

ABERDEEN  PROVING  GROUND,  MD  2100S 
COMMANDER 

US  ARMY  TRAINING  &  DOCTRINE  COMMAND 
ATTN  ATORI-OP-SY1 
PT  MONROE,  VA  23651 

COMMANDER 

WHITE  SANDS  MISSILE  RANGE 
DEPARTMENT  OF  THE  ARMY 
ATTN  STEWS-TE-AN,  J.  OKUMA 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 

OFFICER-IN-CHARGE 
CIVIL  ENGINEERING  LABORATORY 
NAVAL  CONSTRUCTION  BATTALION  CENTER 
ATTN  CODE  L08A  (LIBRARY) 

ATTN  CODE  L08A 

PORT  HUENEME,  CA  93041 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
ATTN  AIR-350F 
WASHINGTON,  DC  21360 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
ATTN  PME  117-215 
WASHINGTON,  DC  20360 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
ATTN  CODE  015,  C.  FLETCHER 
ATTN  CODE  7240,  S.  W.  LICHTMAN 
SAN  DIEGO,  CA  92152 

COMMANDING  OFFICER 
NAVAL  ORDNANCE  STATION 
ATTN  STANDARDIZATION  DIV 
INDIANHEAD,  MD  20640 

SUPERINTENDENT  (CODE  1424) 

NAVAL  POSTGRADUATE  SCHOOL 
ATTN  CODE  1424 
MONTEREY,  CA  93940 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
ATTN  CODE  4104,  EMANUAL  L.  BRANCATO 
ATTN  CODE  2627,  DORIS  R.  FOLEN 
ATTN  CODE  662  3,  RICHARD  L.  STATLER 
ATTN  CODE  6624 
WASHINGTON,  DC  20375 


COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
ATTN  CODE  6174D2,  EDWARD  F.  DUFFY 
WASHINGTON,  DC  20362 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN  CODE  WA51RH,  RM  130-108 
ATTN  CODE  F32,  EDWIN  R.  RATHBURN 
WHITE  OAK,  SILVER  SPRING,  MD  20910 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
ATTN  CODE  DF-56 
DAHLGREN,  VA  22448 

COMMANDER 

NAVAL  WEAPONS  CENTER 
ATTN  CODE  533,  TECH  LIB 
CHINA  LAKE,  CA  93555 

COMMANDING  OFFICER 

NAVAL  WEAPONS  EVALUATION  FACILITY 

KIRTLAND  AIR  FORCE  BASE 

ATTN  CODE  AT- 6 

ALBUQUERQUE,  NM  87117 

OFFICE  OF  NAVAL  RESEARCH 
ATTN  CODE  427 
ARLINGTON,  VA  22217 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 

ATTN  NSP-2701 ,  JOHN  W.  PITSENBERGER 

ATTN  NSP-2342,  RICHARD  L.  COLEMAN 

ATTN  NSP-43,  TECH  LIB 

ATTN  NSP-27334 

ATTN  NSP-230,  D.  GOLD 

WASHINGTON,  DC  20376 

COMMANDER 

AERONAUTICAL  SYSTEMS  DIVISION,  AFSC 
ATTN  ASD-YH-EX 
ATTN  ENFTV 

WRIGHT-PATTERSON  AFB,  OH  45333 

AIR  FORCE  TECHNICAL  APPLICATIONS 
CENTER 

ATTN  TFS,  M.  SCHNEIDER 
PATRICK  AFB,  FL  32925 


73 


DISTRIBUTION  (Cont *d) 


AF  WEAPONS  LABORATORY,  AFSC 
ATTN  NTN 
ATTN  NT 

ATTN  EL,  CARL  E.  BAUM 
ATTN  EL XT 
ATTN  SUL 
ATTN  CA 

ATTN  ELA,  J.  P.  CASTILLO 
ATTN  ELP 

ATTN  ELT,  W,  PAGE 
ATTN  NXS 

KIRTLAND  AFB,  NM  87117 
DIRECTOR 

AIR  UNIVERSITY  LIBRARY 
DEPARTMENT  OF  THE  AIR  FORCE 
ATTN  AUL-LSE-70-250 
MAXWELL  AFB,  AL  36112 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVIS ION/ YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
ATTN  YSEA 

HANSCOM  AFB,  MA  01731 
COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
ATTN  NICD  LIBRARY 
ATTN  ETDP,  B.  L.  BALLARD 
WRIGHT-PATTERSON  AFB,  OH  45433 

COMMANDER 

OGDEN  ALC/MMEDDE 

DEPARTMENT  OF  THE  AIR  FORCE 

ATTN  OO-ALC/MMETH,  P.  W.  BERTHEL 

ATTN  MMEDO,  LEO  KIDMAN 

ATTN  MAJ  R.  BLACKBURN 

HILL  AFB,  UT  84406 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
ATTN  TSLD 

GRIFFISS  AFB,  NY  13441 
COMMANDER 

SACRAMENTO  AIR  LOGISTICS  CENTER 
DEPARTMENT  OF  THE  AIR  FORCE 
ATTN  MMCRS,  H,  A.  PE  LM  ASTRO 
ATTN  MMIRA,  J.  W.  DEMES 
ATTN  MMSREM,  F.  R«  SPEAR 
MCCLELLAN  AFB,  CA  95652 

SAMSO/IN 

AIR  FORCE  SYSTEMS  COMMAND 
PO  BOX  92960 
WORLDWAY  POSTAL  CENTER 
(INTELLIGENCE) 

ATTN  IND 

LOS  ANGELES,  CA  90009 


SAMSO/MN 

AIR  FORCE  SYSTEMS  COMMAND 
(MINUTEMAN) 

ATTN  MNNH,  MAJ  M.  BARAN 

ATTN  MNNH,  CAPT  R.  I.  LAWRENCE 

NORTON  AFB,  CA  92409 

SAMSO/YA 

AIR  FORCE  SYSTEMS  COMMAND 
PO  BOX  92960 
WORLDWAY  POSTAL  CENTER 
ATTN  YAPC 

LOS  ANGELES,  CA  90009 

STRATEGIC  AIR  COMMAND/XPFS 
ATTN  NRI-STINFO  LIBRARY 
ATTN  DEL 

ATTN  GARNET  E.  MATZKE 

ATTN  XPFS,  MAJ  BRIAN  G.  STEPHAN 

OFFUTT  AFB,  NB  68113 

DEPARTMENT  OF  ENERGY 

ALBUQUERQUE  OPERATIONS  OFFICE 

PO  BOX  5400 

ATTN  TECH  LIBRARY 

ATTN  OPERATIONAL  SAFETY  DIV 

ALBUQUERQUE,  NM  87115 

UNIVERSITY  OF  CALIFORNIA 

LAWRENCE  LIVERMORE  LABORATORY 

PO  BOX  808 

ATTN  TECH  INFO  DEPT 

ATTN  L-06,  T.  DONICH 

ATTN  L-545,  D.  MEEKER 

ATTN  L-1 56,  E.  MILLER 

ATTN  L-1 0,  H.  KRUGER 

ATTN  FCPRL 

ATTN  H.  S.  CABAYAN 

LIVERMORE,  CA  94550 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

PO  BOX  1663 

ATTN  BRUCE  W.  NOEL 

ATTN  CLARENCE  BENTON 

LOS  ALAMOS,  NM  87545 

SANDIA  LABORATORIES 
PO  BOX  5800 
ATTN  C.  N.  VITTITOE 
ATTN  R.  L.  PARKER 
ATTN  ELMER  F.  HARTMAN 
ALBUQUERQUE,  NM  87115 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  RM  SG48,  HQ,  BLDG 
( OS I/NED/NWB ) 

WASHINGTON,  DC  20505 


74 


DISTRIBUTION  (Cont  d) 


EfENSE  ELECTRIC  POWER  ADMIN 
EPARTMENT  OF  THE  interior 

nterior  south  bldg,  312 
.TTW  L*  O  *  WEILL 
IASHINGTON,  DC  20240 

DEPARTMENT  of  transportation 

HEADQUARTERS  SEC  DIV,  AS E- 300 
goo  INDEPENDENCE  AVENUE,  SW 
SEC  DIV  ASE-300 


AEROSPACE  CORPORATION 
PO  BOX  92957 
ATTN  C.  B.  PEARLSTON 
ATTN  IRVING  M.  GARPUNKEL 
ATTN  JULIAN  REINHEIMER 
ATTN  LIBRARY 
ATTN  CHARLES  GREENHOW 
LOS  ANGELES,  CA  90009 


AGABIAN  ASSOCIATES 
250  NORTH  NASH  STREET 
ATTN  LIBRARY 
EL  SEGUNDO,  CA  90245 

AVCO  RESEARCH  S  SYSTEMS  GROUP 
201  LOWELL  STREET 
WILMINGTON,  MA  01887 

BATTELLE  MEMORIAL  INSTITUTE 
605  KING  AVENUE 
ATTN  ROBERT  H.  BLAZER 
ATTN  EUGENE  R.  LEACH 
COLUMBUS,  OH  43201 


BOM  CORPORATION 
7915  JONES  BRANCH  DRIVE 
ATTN  CORPORATE  LIBRARY 
MCLEAN,  VA  22102 


BDM  CORPORATION 
PO  BOX  9274 

ALBUQUERQUE  INTERNATIONAL 

ATTN  LIB 

ALBUQUERQUE,  NM  87119 


LEND  IX  CORPORATION 

lESBARCH  LABORATORIES  DIVISIO 

SEND  IX  CENTER 

httn  max  frank 

_ _ _ r*T  n  mt  4B075 


nFNDIX  CORPORATION 
NAVIGATION  AND  CONTROL  GROUP 


ATTN  DEPT  6401 
teterboro,  nj 


07608 


BOEING  COMPANY 
PO  BOX  3707 

ATTN  HOWARD  W.  WICKLEIN 
ATTN  D.  E.  ISBELL 
ATTN  DAVID  KEMLE 
ATTN  B.  C.  HANRAHAN 
ATTN  KENT  TECH  LIB 
SEATTLE,  WA  98124 

BROWN  ENGINEERING  COMPANY,  INC 
CUMMINGS  RESEARCH  PARK 
ATTN  FRED  LEONARD 
HUNTSVILLE,  AL  35807 


LURROUGHS  CORPORATION 
•EDERAL  AND  SPECIAL  SYSTEMS  GROUP 
CENTRAL  AVE  AND  ROUTE  252 
PO  BOX  517 


CALSPAN  CORPORATION 
PO  BOX  400 
ATTN  TECH  LIBRARY 
BUFFALO,  NY  14225 

CHARLES  STARK  DRAPER  LABORATORY,  INC 

555  TECHNOLOGY  SQUARE 

ATTN  KENNETH  FERTIG 

ATTN  TIC  MS  74 

CAMBRIDGE,  MA  021 39 


CINCINNATI  ELECTRONICS  CORPORATION 
2630  GLENDALE  -  MILFORD  ROAD 


COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVD 
ATTN  RAMONA  BRIGGS 
FALLS  CHURCH,  VA  22046 


COMPUTER  SCIENCES  CORPORATION 
1400  SAN  MATEO  BLVD,  SE 
ATTN  RICHARD  H.  DICKHAUT 
ATTN  ALVIN  SCHIFF 
ALBUQUERQUE,  NM  87108 


CONTROL  DATA  CORPORATION 
PO  BOX  0 

ATTN  JACK  MEEHAN 
MINNEAPOLIS,  MN  55440 


CUTLER-HAMMER,  INC 
AIL  DIVISION 
COMAC  ROAD 
ATTN  EDWARD  KARPEN 
DEER  PARK,  NY  1 1 729 


75 


DISTRIBUTION  (Cont'd) 


DIKEWOOD  CORPORATION 
1613  UNIVERSITY  BLVD,  NE 
ATTN  TECH  LIB 
ATTN  L.  WAYNE  DAVIS 
ALBUQUERQUE,  NM  87102 

DIKEWOOD  CORPORATION 

2716  OCEAN  &  PARK  BLVD,  SUITE  3000 

ATTN  K.  LEE 

SANTA  MONICA,  CA  90405 

E-SYSTEMS,  INC 
GREENVILLE  DIVISION 
PO  BOX  1056 
ATTN  JOLETA  MOORE 
GREENVILLE,  TX  75401 

EFFECTS  TECHNOLOGY,  INC 
5383  HOLLISTER  AVENUE 
ATTN  S.  CLOW 
SANTA  BARBARA,  CA  93111 

EG&G  WASHINGTON  ANALYTICAL  SERVICES 
CENTER,  INC 
PO  BOX  10218 
ATTN  C.  GILES 
ALBUQUERQUE,  NM  87114 

EXXON  NUCLEAR  COMPANY,  INC 
RESEARCH  AND  TECHNOLOGY  CENTER 
2955  GEORGE  WASHINGTON  WAY 
ATTN  DR.  A.  W.  TRIVELPIECE 
RICHLAND,  WA  99352 

FAIRCHILD  CAMERA  &  INSTRUMENT  CORP 
464  ELLIS  STREET 
ATTN  DAVID  K.  MYERS 
MOUNTAIN  VIEW,  CA  94040 

FORD  AEROSPACE  &  COMMUNICATIONS  CORP 

3939  FABIAN  WAY 

ATTN  TECHNICAL  LIBRARY 

PALO  ALTO,  CA  94303 

FORD  AEROSPACE  &  COMMUNICATIONS 
OPERATIONS 

FORD  &  JAMBOREE  ROADS 
ATTN  KEN  C.  ATTINGER 
NEWPORT  BEACH,  CA  92663 

FRANKLIN  INSTITUTE 
20TH  STREET  &  PARKWAY 
ATTN  RAMIE  H.  THOMPSON 
PHILADELPHIA,  PA  19103 

GENERAL  DYNAMICS  CORP 
ELECTRONICS  DIVISION 
PO  BOX  81 1 25 
ATTN  RSCH  LIB 
SAN  DIEGO,  CA  92138 


GENERAL  DYNAMICS  CORPORATION 

INTER-DIVISION  RESEARCH  LIBRARY 

KEARNY  MESA 

PO  BOX  80847 

ATTN  RSCH  LIB 

SAN  DIEGO,  CA  92138 

GENERAL  ELECTRIC  CO. -TEMPO 
CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (PO  DRAWER  QQ) 
ATTN  DASIAC 

ATTN  ROYDEN  R.  RUTHERFORD 
ATTN  WILLIAM  MC NAM ERA 
SANTA  BARBARA,  CA  93102 

GENERAL  ELECTRIC  COMPANY 
AEROSPACE  ELECTRONICS  SYSTEMS 
FRENCH  ROAD 

ATTN  CHARLES  M.  HEWISON 
UTICA,  NY  73503 

GENERAL  ELECTRIC  COMPANY 
PO  BOX  5000 
ATTN  TECH  LIB 
BINGHAMTON,  NY  1 3902 

GENERAL  ELECTRIC  CO. -TEMPO 
ALEXANDRIA  OFFICE 
HUNTINGTON  BUILDING,  SUITE  300 
2560  HUNTINGTON  AVENUE 
ATTN  DASIAC 
ALEXANDRIA,  VA  22303 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
PO  BOX  6770 
ATTN  TECH  INFO  OFFICE 
SANTA  BARBARA,  CA  93111 

GEORGIA  INSTITUTE  OF  TECHNOLOGY 
GEORGIA  TECH  RESEARCH  INSTITUTE 
ATTN  R.  CURRY 
ATLANTA,  GA  30332 

GEORGIA  INSTITUTE  OF  TECHNOLOGY 
OFFICE  OF  CONTRACT  ADMINISTRATION 
ATTN  RSCH  SECURITY  COORDINATOR 
FOR  HUGH  DENNY 
ATLANTA,  GA  30332 

GRUMMAN  AEROSPACE  CORPORATION 
SOUTH  OYSTER  BAY  ROAD 
ATTN  L-01  35 
BETHPAGE,  NY  11714 


76 


DISTRIBUTION  (Cont'd) 


GTE  SYLVAN IA,  INC 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 

ATTN  CHARLES  A.  THORNHILL,  LIBRARIAN 
ATTN  LEONARD  L.  BLAISDELL 
NEEDHAM  HEIGHTS,  MA  02194 

GTE  SYLVAN IA,  INC 
189  B  STREET 

ATTN  CHARLES  H.  RAMS BOTTOM 

ATTN  DAVID  D.  FLOOD 

ATTN  EMIL  P.  MOTCHOK 

ATTN  H  &  V  GROUP,  MARIO  A.  NUREFORA 

ATTN  J.  WALDRON 

NEEDHAM  HEIGHTS,  MA  02194 

HARRIS  CORPORATION 

HARRIS  SEMICONDUCTOR  DIVISION 

PO  BOX  883 

ATTN  V  PRES  &  MGR  PRGMS  DIV 
MELBOURNE,  FL  32901 

HAZELTINE  CORPORATION 
PULASKI  ROAD 

ATTN  TECH  INFO  CTR,  M.  WAITE 
GREENLAWN,  NY  11740 

HONEYWELL  INC 
AVIONICS  DIVISION 
2600  RIDGEWAY  PARKWAY 
ATTN  S&RC  LIB 
ATTN  RONALD  R.  JOHNSON 
MINNEAPOLIS,  MN  55413 

HONEYWELL  INC 

AVIONICS  DIVISION 

13350  U.S.  HIGHWAY  19  NORTH 

ATTN  MS  725-5,  STACEY  H.  GRAFF 

ATTN  W.  E.  STEWART 

ST  PETERSBURG,  FL  33733 

HUGHES  AIRCRAFT  COMPANY 
CENTINELA  &  TEALE 
ATTN  JOHN  B.  SINGLETARY 
ATTN  CTDC  6/El  1 0 
ATTN  KENNETH  R.  WALKER 
CULVER  CITY,  CA  90230 

I IT  RESEARCH  INSTITUTE 
ELECTROMAG  COMPATABILITY  ANAL  CTR 
NORTH  SEVERN 
ATTN  ACOAT 
ANNAPOLIS,  MD  21402 

I IT  RESEARCH  INSTITUTE 
10  WEST  35TH  STREET 
ATTN  IRVING  N.  MINDEL 
ATTN  JACK  E.  BRIDGES 
CHICAGO,  IL  60616 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ATTN  TECH  INFO  SERVICES 
ARLINGTON,  VA  22202 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
ATTN  TECHNICAL  LIBRARY 
ATTN  ALEXANDER  T.  RICHARDSON 
NUTLEY,  NJ  07110 

IRT  CORPORATION 
PO  BOX  81087 
ATTN  C.  B.  WILLIAMS 
ATTN  DENNIS  SWIFT 
SAN  DIEGO,  CA  92138 

JAYCOR 

11011  TORRE YANA  ROAD 
PO  BOX  85154 
ATTN  ERIC  P.  WENAAS 
ATTN  RALPH  H.  STAHL 
SAN  DIEGO,  CA  92138 

JAYCOR 

205  S.  WHITTING  STREET,  SUITE  500 
ATTN  LIB 

ALEXANDRIA,  VA  22304 

KAMAN  SCIENCES  CORPORATION 

1500  GARDEN  OF  THE  GODS  ROAD 

ATTN  ALBERT  P.  BRIDGES 

ATTN  W.  FOSTER  RICH 

ATTN  WALTER  E.  WARE 

ATTN  FRANK  H.  SHELTON 

ATTN  JERRY  I.  LUBELL 

ATTN  PHIL  TRACY 

COLORADO  SPRINGS,  CO  80907 

LITTON  SYSTEMS,  INC 
DATA  SYSTEMS  DIVISION 
8000  WOODLEY  AVENUE 
ATTN  EMC  GP 
ATTN  M848-61 
VAN  NUYS,  CA  91409 

LITTON  SYSTEMS,  INC 
AMERCOM  DIVISION 
5115  CALVERT  ROAD 
ATTN  J.  SKAGGS 
COLLEGE  PARK,  MD  20740 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 

PO  BOX  504 

ATTN  I.  ROSSI 

ATTN  SAMUEL  I.  TAIMUTY 

ATTN  H.  E.  THAYN 

ATTN  GEORGE  F.  HEATH 

ATTN  BENJAMIN  T.  KIMURA 

SUNNYVALE,  CA  94086 


77 


DISTRIBUTION  (Cont'd) 


) 

* 


LOCKHEED  MISSILE  &  SPACE  CO,  INC 
3251  HANOVER  STREET 
ATTN  TECH  INFO  CTR  D/COLL 
PALO  ALTO,  CA  94304 

MIT  LINCOLN  LABORATORY 
PO  BOX  73 

ATTN  LEONA  LOUGHLIN 
LEXINGTON,  MA  02173 

MARTIN  MARIETTA  CORPORATION 
ORLANDO  DIVISION 
PO  BOX  5837 
ATTN  MONA  C.  GRIFFITH 
ORLANDO,  FL  32805 

MCDONNELL  DOUGLAS  CORPORATION 
PO  BOX  516 
ATTN  TOM  ENDER 
ST  LOUIS,  MO  63166 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
ATTN  STANLEY  SCHNEIDER 
ATTN  TECH  LIBRARY  SERVICES 
HUNTINGTON  BEACH,  CA  92647 

MISSION  RESEARCH  CORPORATION 

PO  DRAWER  719 

ATTN  EMP  GROUP 

ATTN  WILLIAM  C.  HART 

ATTN  C.  LONGMIRE 

SANTA  BARBARA,  CA  93102 

MISSION  RESEARCH  CORPORATION 
PO  BOX  7816 
ATTN  WERNER  STARK 
ATTN  ROY  STRAYER,  JR 
COLORADO  SPRINGS,  CO  80933 

MISSION  RESEARCH  CORPORATION 
EM  SYSTEM  APPLICATIONS  DIVISION 
1400  SAN  MATEO  BLVD,  SE,  SUITE  A 
ATTN  DAVID  E.  MEREWETHER 
ALBUQUERQUE,  NM  87108 

MISSION  RESEARCH  CORPORATION- 
SAN  DIEGO 
PO  BOX  1209 

ATTN  V.  A.  J.  VAN  LINT 
LA  JOLLA,  CA  92038 

MITRE  CORPORATION 
PO  BOX  208 

ATTN  M.  F.  FITZGERALD 
BEDFORD,  MA  01730 


NORDEN  SYSTEMS,  INC 
HELEN  STREET 
ATTN  TECHNICAL  LIBRARY 
NORWALK,  CT  06856 

NORTHROP  RESEARCH  TECHNOLOGY  CENTER 
ONE  RESEARCH  PARK 
ATTN  LIBRARY 

PALOS  VERDES  PENN,  CA  90274 

NORTHROP  CORPORATION 
ELECTRONIC  DIVISION 
2301  WEST  1 20TH  STREET 
ATTN  LEW  SMITH 
ATTN  RAD  EFFECTS  GRP 
HAWTHORNE,  CA  90250 

PHYSICS  INTERNATIONAL  COMPANY 
2700  MERCED  STREET 
ATTN  DOC  CON 
SAN  LEANDRO,  CA  94577 

R&D  ASSOCIATES 

PO  BOX  9695 

ATTN  S.  CLAY  ROGERS 

ATTN  CHARLES  MO 

ATTN  RICHARD  R.  SCHAEFER 

ATTN  DOC  CON 

ATTN  M.  GROVER 

ATTN  C.  MACDONALD 

ATTN  J.  BOMBARDT 

MARINA  DEL  REY,  CA  90291 

RAND  CORPORATION 

1700  MAIN  STREET 

ATTN  LIB-D 

ATTN  W.  SOLLFREY 

SANTA  MONICA,  CA  90406 

RAYTHEON  COMPANY 
HARTWELL  ROAD 
ATTN  GAJANAN  H,  JOSHI 
BEDFORD,  MA  01730 

RAYTHEON  COMPANY 
528  BOSTON  POST  ROAD 
ATTN  HAROLD  L.  FLESCHER 
SUDBURY,  MA  01776 

RCA  CORPORATION 
GOVERNMENT  SYSTEMS  DIVISION 
ASTRO  ELECTRONICS 
PO  BOX  80 0,  LOCUST  CORNER 
EAST  WINDSOR  TOWNSHIP 
ATTN  GEORGE  J.  BRUCKER 
PRINCETON,  NJ  08540 


78 


DISTRIBUTION  (Cont'd) 


RCA  CORPORATION 

DAVID  SARNOFF  RESEARCH  CENTER 

PO  BOX  432 

ATTN  SECURITY  DEPT,  L.  MINICH 
PRINCETON,  NJ  08540 

RCA  CORPORATION 
CAMDEN  COMPLEX 
FRONT  &  COOPER  STREETS 
ATTN  OLIVE  WHITEHEAD 
ATTN  R.  W.  ROSTROM 
CAMDEN,  NJ  08012 

ROCKWELL  INTERNATIONAL  CORP 

PO  BOX  3105 

ATTN  N.  J.  RUDIE 

ATTN  J.  L.  MONROE 

ATTN  V.  J.  MICHEL 

ATTN  D/243-068,  031-CA31 

ANAHEIM,  CA  92803 

ROCKWELL  INTERNATIONAL  CORP 
SPACE  DIVISION 

12214  SOUTH  LAKEWOOD  BOULEVARD 
ATTN  B.  E.  WHITE 
DOWNEY,  CA  90241 

ROCKWELL  INTERNATIONAL  CORPORATION 

81 5  LA  PHAM  STREET 

ATTN  B-1,  DIV  TIC  ( BAOB ) 

EL  SEGUNDO,  CA  90245 

ROCKWELL  INTERNATIONAL  CORPORATION 
PO  BOX  369 
ATTN  F.  A.  SHAW 
CLEARFIELD,  UT  8401 5 

SANDERS  ASSOCIATES,  INC 
95  CANAL  STREET 

ATTN  1-6270,  R.  G.  DESPATHY,  SR  P  E 
NASHUA,  NH  03060 

SCIENCE  APPLICATIONS,  INC 
PO  BOX  277 

ATTN  FREDERICK  M.  TESCHE 
BERKELEY,  CA  94701 

SCIENCE  APPLICATIONS,  INC 
PO  BOX  2351 
ATTN  R.  PARKINSON 
LA  JOLLA,  CA  92038 

SCIENCE  APPLICATIONS,  INC 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 
ATTN  NOEL  R.  BYRN 
HUNTSVILLE,  AL  35805 


SCIENCE  APPLICATIONS,  INC 
8400  WESTPARK  DRIVE 
ATTN  WILLIAM  L.  CHADSEY 
MCLEAN,  VA  22101 

SINGER  COMPANY 
1150  MC  BRIDE  AVENUE 
ATTN  TECH  INFO  CTR 
LITTLE  FALLS,  NJ  07424 

SPERRY  RAND  CORPORATION 
SPERRY  MICROWAVE  ELECTRONICS 
PO  BOX  4648 
ATTN  MARGARET  CORT 
CLEARWATER,  FL  33518 

SPERRY  RAND  CORPORATION 
SPERRY  DIVISION 
MARCUS  AVENUE 
ATTN  TECH  LIB 
GREAT  NECK,  NY  11020 

SPERRY  RAND  CORPORATION 
SPERRY  FLIGHT  SYSTEMS 
PO  BOX  21111 
ATTN  D.  ANDREW  SCHOW 
PHOENIX,  AZ  85036 

SPIRE  CORPORATION 
PO  BOX  D 

ATTN  JOHN  R.  UGLUM 
ATTN  ROGER  G.  LITTLE 
BEDFORD,  MA  01730 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
ATTN  ARTHUR  LEE  WHITSON 
MENLO  PARK,  CA  94025 

SYSTEMS,  SCIENCE  &  SOFTWARE,  INC 
PO  BOX  1620 
ATTN  ANDREW  R.  WILSON 
LA  JOLLA,  CA  92038 

TEXAS  INSTRUMENTS,  INC 
PO  BOX  601 5 
ATTN  TECH  LIB 
ATTN  DONALD  J  MANUS 
DALLAS,  TX  75265 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

ATTN  0.  E.  ADAMS 

ATTN  R.  K.  PLEBUCH 

ATTN  L.  R.  MAGNOLIA 

ATTN  H.  H,  HOLLOWAY 

ATTN  W.  GARGARO 

REDONDO  BEACH,  CA  90278 


79 


•  dr> 


DISTRIBUTION  (Cont'd) 


TEXAS  TECH  UNIVERSITY 
PO  BOX  5404  NORTH  COLLEGE  STATION 
ATTN  TRAVIS  L.  SIM»S0N 
LUBBOCK,  TX  79417 

UNITED  TECHNOLOGIES  CGRP 
HAMILTON  STANDARD  DIVISION 
BRADLEY  INTERNATIONAL  AIRPORT 
ATTN  CHIEF  ELEC  DESIGN 
WINDSOR  LOCKS,  CT  06069 

WESTINGHOUSE  ELECTRIC  CORPORATION 
ADVANCED  ENERGY  SYSTEMS  DIV 
PO  BOX  10864 
ATTN  TECH  LIB 
PITTSBURGH,  PA  15236 

MDB  ENTERPRISES  INC 
9101  BARTON 

ATTN  MONTI  R.  WILSON  (10  COPIES) 

OVERLAND  PARK,  KS  6621 4 

US  ARMY  ELECTRONICS  RESEARCH 
&  DEVELOPMENT  COMMAND 
ATTN  TECHNICAL  DIRECTOR,  DRDEL-CT 
ATTN  R.  HARMAN,  DRDEL-MA 

HARRY  DIAMOND  LABORATORIES 

ATTN  CO/TD/TSO/D IVISION  DIRECTORS 

ATTN  RECORD  COPY,  81200 

ATTN  HOL  LIBRARY,  81100  (2  COPIES) 

ATTN  HOL  LIBRARY,  81100  (WOODBRIDGE) 

ATTN  TECHNICAL  REPORTS  BRANCH,  81300  (3  COPIES) 

ATTN  LEGAL  OFFICE,  97000 

ATTN  CHAIRMAN,  EDITORIAL  COMMITTEE 

ATTN  CHIEF,  13000 

ATTN  CHIEF,  21000 

ATTN  CHIEF,  21100 

ATTN  CHIEF,  21200 

ATTN  CHIEF,  21300  (10  COPIES) 

ATTN  CHIEF,  21400 
ATTN  CHIEF,  21500 
ATTN  CHIEF,  22000 
ATTN  CHIEF,  22100 
ATTN  CHIEF,  22300 
ATTN  CHIEF,  22800 
ATTN  CHIEF,  22900 
ATTN  CHIEF,  20240 
ATTN  HUTTLIN,  G. ,  22900 
ATTN  CHASE,  R,  ,  21100 
ATTN  VALLIN,  J. ,  22100 
ATTN  FEMENIAS,  R,  ,  22100 
ATTN  LEPOER,  K.  ,  21100 
ATTN  REYZER,  R  ,  21300 
ATTN  LANHAM,  C. ,  00213 
ATTN  PEPERONE,  S. ,  36000 
ATTN  RESSLER,  M.  ,  15300 
ATTN  LEWIS,  D. ,  15300 
ATTN  SHREVE,  J. ,  15300 
ATTN  BIEHL,  B.  ,  15300 
ATTN  SEBOL,  E.  ,  36200 
ATTN  SANN,  X.,  15000 
ATTN  IMMERMAN,  P.  ,  1 5400 
ATTN  MEYER,  O. ,  22800 
ATTN  SOLN,  J. ,  22300 
ATTN  CUNEO,  A.,  21400 


80 


